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Technology Focus: Composites & Ceramics 


$ Improving Thermomechanical Properties of SiC/SiC Composites 

A heat treatment increases thermal conductivity and increases creep resistance. 

John H. Glenn Research Center, Cleveland, Ohio 
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Figure 1. Room-Temperature Thermal Conductivities were measured on unannealed and annealed 
specimens of composite panels containing SiC matrices formed in various ways as described in the last 
paragraph of the main text. 
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Temperature of 1,315 °C and 
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Figure 2. Creep Strains were measured on unannealed and annealed specimens of composite panels 
containing SiC matrices like two of those represented in Figure 1. 


Today, a major thrust toward improv- 
ing the thermomechanical properties of 
engine components lies in the develop- 
ment of fiber-reinforced silicon carbide 
matrix composite materials, including 
SiC-fiber/SiC-matrix composites. These 
materials are lighter in weight and capa- 
ble of withstanding higher tempera- 
tures, relative to state-of-the-art metallic 
alloys and oxide-matrix composites for 
which maximum use temperatures are 
in the vicinity of 1,100 °C. In addition, 
the toughness or damage tolerance of 
the SiC-matrix composites is signifi- 
cantly greater than that of unreinforced 
silicon-based monolithic ceramics. 

For successful application in advanced 
engine systems, the SiC-matrix compos- 
ites should be able to withstand compo- 
nent service stresses and temperatures 
for the desired component lifetimes. 
Inasmuch as the high-temperature struc- 
tural lives of ceramic materials are typi- 
cally limited by creep-induced growth of 
flaws, a key property required of such 
composite materials is high resistance to 
creep under conditions of use. Also, the 
thermal conductivity of the materials 
should be as high as possible so as to 
minimize component thermal gradients 
and thermal stresses. 

A state-of-the-art SiC-matrix composite 
is typically fabricated in a three-step 
process: (1) fabrication of a component- 
shaped architectural preform reinforced 
by thermally stable high-performance 
fibers, (2) chemical-vapor infiltration 
(CVI) of a fiber-coating material such as 
boron nitride (BN) into the preform, and 
(3) infiltration of an SiC-based matrix into 
the remaining porosity in the preform. 
Generally, the matrices of the highest-per- 
forming composites are fabricated by ini- 
tial use of a CVI SiC matrix component 
that is typically more thermally stable and 
denser than matrix components formed 
by processes other than CVI. As such, the 
initial SiC matrix component made by CVI 
provides better environmental protection 
to the coated fibers embedded within it. 
Also, the denser CVI SiC imparts to the 
composite better resistance to propaga- 
tion of cracks, enhanced thermal conduc- 


tivity, and higher creep resistance. 

For adequate deposition of CVI SiC 
matrix material into a preform, it is typi- 
cally necessary to perform the infiltra- 
tion at a temperature below 1,100 °C. 
The use of an infiltration temperature in 
this range typically results in the forma- 
tion of an SiC matrix that is fairly dense, 
but contains metastable atomic defects 
and is non-stoichiometric in that it con- 
tains a small excess of silicon. These de- 
fects typically exist at the matrix grain 


boundaries, where they can act as scat- 
tered for thermal phonons and degrade 
matrix creep resistance by enhancing 
grain-boundary sliding. 

In order to make it possible to elimi- 
nate these defects and improve the ther- 
momechanical properties of ceramic 
composites that contain SiC matrices, re- 
searchers at Glenn Research Center 
(GRC) have developed a high-tempera- 
ture heat or annealing treatment that 
can be performed after deposition of a 
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CVI SiC matrix into a fiber preform. 
Using (1) SiC fibers of a type developed 
by GRC and denoted “Sylramic-iBN” and 
(2) BN-based fiber coatings, which are 
both stable in their functions under the 
treatment conditions, the GRC re- 
searchers have observed minimal loss of 
strength in composite panels formed 
from two-dimensional architectural pre- 
forms and various contents of CVI SiC. 
More importantly, the NASA treatment 
significantly increased panel thermal 
conductivity and creep resistance, as in- 
dicated in Figures 1 and 2, respectively. 

For the treated panels, the volume 
fractions of the fibers and BN coatings 
were about 36 percent and 8 percent, 
respectively. The first case represented 
in Figure 1 is that of a panel containing 
50 volume percent SiC formed by CVI 
only. In the second case, the remaining 
open porosity in the 35-percent CVI SiC 
matrix was filled by a process denoted 


PIP, which involved repeated infiltration 
and pyrolysis of an SiC-yielding polymer 
(hybrid CVI + PIP). For the third case 
represented in Figure 1, the remaining 
open porosity in a 35-percent CVI SiC 
matrix was filled by a process that in- 
volved repeated melt infiltration (MI) 
of silicon at a temperature near 1,400 
°C (hybrid CVI + MI) . For the first two 
cases, the heat treatment was performed 
after final matrix formation; for the last 
case, the treatment was performed after 
CVI and before MI. In all cases, panel 
thermal conductivity was significantly 
increased by the NASA heat treatment. 
On an absolute scale, Figure 1 also 
shows the detrimental effect of trapped 
porosity on composite conductivity for 
the matrices formed by CVI only or by 
hybrid CVI + PIP, and the beneficial ef- 
fects of the hybrid CVI + MI approach 
that more effectively fills the CVI SiC 
pores. However, Figure 2 shows that in 


the case of the matrix fully formed by 
CVI only, creep resistance and tempera- 
ture capability were greater than in the 
case of the matrix formed by hybrid CVI 
+ MI. Thus, because of their lack of ele- 
mental silicon, the long-term use tem- 
perature of the CVI-only matrix (as well 
as the hybrid CVI + PIP matrix) could 
exceed 1,400 °C, which is well above the 
long-term use temperature for panels 
containing matrices made by hybrid 
CVI + MI. 

This work was done by James A. DiCarlo of 
Glenn Research Center and Ramakrishna 
T. Bhatt of the Army Research Laboratory. 
Further information is contained in a TSP 
(see page 1 ). 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steve Fedor, Mail 
Stop 4-8, 21 000 Brookpark Road, Cleveland, 
Ohio 44135. Refer to LEW-17595-1. 


Aerogel/Particle Composites for Thermoelectric Devices 

Shrinkage is reduced through addition of titania powder. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 



Shrinkages of Representative Aerogels of various densities and of aerogel/titania composites were 
measured. 


Optimizing solution chemistry and 
the addition of titania and fumed silica 
powder reduces shrinkage. These mate- 
rials would serve to increase thermal ef- 
ficiency by providing thermal insulation 
to suppress lateral heat leaks. They 
would also serve to prolong operational 
lifetime by suppressing sublimation of 
certain constituents of thermoelectric 
materials (e.g., sublimation of Sb from 
CoSb 3 ) at typical high operating temper- 
atures. [The use of pure silica aerogels 
as cast-in-place thermal-insulation and 
sublimation-suppression materials was 
described in “Aerogels for Thermal In- 
sulation of Thermoelectric Devices” 
(NPO-40630), NASA Tech Briefs, Vol. 30, 
No. 7 (July 2006), page 50.] 

A silica aerogel is synthesized in a sol- 
gel process that includes preparation of 
a silica sol, gelation of the sol, and dry- 
ing of the gel in a solvent at a supercriti- 
cal temperature and pressure. The util- 
ity of pure silica aerogel is diminished by 
a tendency to shrink (and, therefore, 
also to crack) during the gelation and 
supercritical-drying stages. Moreover, to 
increase suppression of sublimation, it is 
advantageous to make an aerogel having 
greater density, but shrinkage and crack- 
ing tend to increase with density. 

A composite material of the type 


under investigation consists mostly of ti- 
tania oxide powder particles and a small 
addition of fumed silica powder, which 
are mixed into the sol along with other 


ingredients prior to the gelation stage of 
processing. The silica aerogel and 
fumed silica act as a binder, gluing the ti- 
tania particles together. It is believed 
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that the addition of fumed silica stiffens 
the aerogel network and reduces shrink- 
age during the supercritical-drying 
stage. Minimization of shrinkage en- 
ables establishment of intimate contact 
between thermoelectric legs and the 
composite material, thereby maximizing 
the effectiveness of the material for ther- 
mal insulation and suppression of subli- 
mation. 

To some extent, the properties of the 
composite can be tailored via the propor- 
tions of titania and other ingredients. In 
particular (see figure), the addition of a 


suitably large proportion of titania (e.g., 
0.6 g/ cm 3 ) along with a 10-percent in- 
crease in the amount of tetraethy- 
lorthosilicate [TEOS (an ingredient of 
the sol)] to an aerogel component hav- 
ing a density 40 mg/cm 3 makes it possible 
to cast a high-average-density (>0.1 
g/ cm 3 ) aerogel/particle composite hav- 
ing low shrinkage (2.3 percent) . 

This work was done by Jong-Ah Paik, Jef- 
frey Sakamoto, and Steven Jones of Caltech 
for NASA’s Jet Propulsion Laboratory. 
Further information is contained in a TSP 
(see page 1 ). 


In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Innovative Technology Assets Management 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-42031, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 


Patches for Repairing Ceramics and Ceramic-Matrix Composites 

Patches are simply pressed in place, then heated. 

Lyndon B. Johnson Space Center, Houston, Texas 


Patches consisting mostly of ceramic 
fabrics impregnated with partially cured 
polymers and ceramic particles are 
being developed as means of repairing 
ceramics and ceramic-matrix compos- 
ites (CMCs) that must withstand tem- 
peratures above the melting points of 
refractory metal alloys. These patches 
were conceived for use by space-suited, 
space-walking astronauts in repairing 
damaged space-shuttle leading edges: as 
such, these patches could be applied in 
the field, in relatively simple proce- 
dures, and with minimal requirements 
for specialized tools. These design char- 
acteristics also make the patches useful 
for repairing ceramics and CMCs in ter- 
restrial settings. 

In a typical patch as supplied to an as- 
tronaut or repair technician, the poly- 
mer would be in a tacky condition, de- 
noted as an “A” stage, produced by 
partial polymerization of a monomeric 
liquid. The patch would be pressed 


against the ceramic or CMC object to be 
repaired, relying on the tackiness for 
temporary adhesion. The patch would 
then be bonded to the workpiece and 
cured by using a portable device to heat 
the polymer to a curing temperature 
above ambient temperature but well 
below the maximum operating tempera- 
ture to which the workpiece is expected 
to be exposed. The patch would subse- 
quently become pyrolized to a 
ceramic/glass condition upon initial ex- 
posure to the high operating tempera- 
ture. In the original space-shuttle appli- 
cation, this exposure would be 
Earth-atmosphere-reentry heating to 
about 3,000 °F (about 1,600 °C). 

Patch formulations for space-shuttle 
applications include SiC and Zr0 2 fab- 
rics, a commercial SiC-based pre-ceramic 
polymer, and suitable proportions of 
both SiC and Zr0 2 particles having sizes 
of the order of 1 pm. These formulations 
have been tailored for the space-shuttle 


leading-edge material, atmospheric com- 
position, and reentry temperature profile 
so as to enable repairs to survive re-entry 
heating with expected margin. Other for- 
mulations could be tailored for specific 
terrestrial applications. 

This work was done by Peter A. Hogenson, 
Gordon R. Toombs, Steven Adam, and James 
V. Tompkins of The Boeing Co. for Johnson 
Space Center. 

Title to this invention has been waived 
under the provisions of the National Aero- 
nautics and Space Act {42 U.S.C. 2457(f)}, 
to The Boeing Company. Inquiries concern- 
ing licenses for its commercial development 
should be addressed to: 

The Boeing Company 
PO Box 2515 
2201 Seal Beach Blvd. 

Seal Beach, CA 90740-1515 
Phone No. (562) 797-2020 
Refer to MSC-2401 8, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 


Lower-Conductivity Ceramic Materials for Thermal-Barrier 
Coatings 

Thermal conductivities of certain pyrochlore oxides can be reduced by doping. 

John H. Glenn Research Center, Cleveland, Ohio 


Doped pyrochlore oxides of a type de- 
scribed below are under consideration as al- 
ternative materials for high-temperature 
thermal-barrier coatings (TBCs). In com- 
parison with partially-yttria-stabilized zirco- 
nia (YSZ), which is the state-of-the-art TBC 


material now in commercial use, these 
doped pyrochlore oxides exhibit lower ther- 
mal conductivities, which could be ex- 
ploited to obtain the following advantages: 

• For a given difference in temperature 
between an outer coating surface and 


the coating/ substrate interface, the 
coating could be thinner. Reductions in 
coating thicknesses could translate to re- 
ductions in weight of hot-section com- 
ponents of turbine engines (e.g., com- 
bustor liners, blades, and vanes) to 
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Measured Thermal Conductivities of hot-pressed 
tions illustrate the benefit afforded by doping. 

which TBCs are typically applied. 

• For a given coating thickness, the differ- 
ence in temperature between the outer 
coating surface and the coating/sub- 
strate interface could be greater. For tur- 
bine engines, this could translate to 
higher operating temperatures, with 
consequent increases in efficiency and 
reductions in polluting emissions. 

TBCs are needed because the temper- 
atures in some turbine-engine hot sec- 
tions exceed the maximum tempera- 
tures that the substrate materials 
(superalloys, Si-based ceramics, and oth- 
ers) can withstand. YSZ TBCs are ap- 
plied to engine components as thin lay- 
ers by plasma spraying or electron-beam 
physical vapor deposition. During oper- 
ation at higher temperatures, YSZ layers 


disks of several related pyrochlore oxide composi- 

undergo sintering, which increases their 
thermal conductivities and thereby ren- 
ders them less effective as TBCs. More- 
over, the sintered YSZ TBCs are less tol- 
erant of stress and strain and, hence, are 
less durable. 

The materials that are sought as alter- 
natives to YSZ are required to have and 
retain lower thermal conductivities and 
to be better able to withstand tempera- 
tures that degrade TBCs made of YSZ. 
The undoped versions of the doped py- 
rochlore oxides of the type now under 
consideration as alternatives to YSZ are 
of general composition Ma 2 Mb 2 O7, 
where Ma denotes a 3+ cation (for exam- 
ple, La to Lu) and Mb a 4+ cation (for 
example, Zr, Hf, Ti). Doping has been 
investigated as a means of reducing ther- 


mal conductivities even further below 
those of YSZ coatings. In the doping ap- 
proach investigated thus far, another 
cation is substituted for part of Ma, yield- 
ing a general composition of Ma^_ x 
M x Mb 2 O7, where xlies between 0 and 0.5 
and M denotes a rare-earth or other suit- 
able element. 

In experiments, powders of various 
compositions were synthesized by a 
modified sol-gel method and calcined at 
appropriate temperatures to convert 
them into compounds of pyrochlore 
structure as confirmed by x-ray diffrac- 
tion. These powders were hot pressed 
into dense disks of 1-in. (2.54-cm) diam- 
eter. The thermal conductivities of the 
disks were measured at various tempera- 
tures up to 1,550 °C by use of a steady- 
state laser heat-flux technique. The fig- 
ure presents results of such 
measurements performed on several 
materials of general composition 
La 2 _ x (Gd and/or Yb) x Zr 2 0 7 , where x = 0 
or 0.3. The thermal conductivities of all 
doped samples (x = 0.3) were less than 
those of the undoped (x = 0) sample 
[La 2 Zr 2 C>7]. The lowest conductivity — 
ranging from 40 to 50 percent below 
that of undoped sample — was exhib- 
ited by the sample co-doped with both 
Gd andYb. 

This work was done by Narottam P. Bansal 
of Glenn Research Center and Dongming 
Zhu of the U. S. Army Research Laboratory. 
Further information is contained in a TSP (see 
page 1). 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed to: 
NASA Glenn Research Center 
Innovative Partnerships Office 
Attn: Steve Fedor 
Mail Stop 4-8 
21000 Brookpark Road 
Cleveland, Ohio 44135 
Refer to LEW-17469-1. 
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Electronics/Computers 


An Alternative for Emergency Preemption of Traffic Lights 

This system resolves potential conflicts among emergency vehicles. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


An electronic communication-and- 
control system has been developed as a 
prototype of advanced means of auto- 
matically modifying the switching of 
traffic lights to give priority to emer- 
gency vehicles. This system could be 
used alternatively or in addition to 
other emergency traffic-light-preemp- 
tion systems, including a variety of sys- 
tems now in use as well as two proposed 
systems described in “Systems Would 
Preempt Traffic Lights for Emergency 
Vehicles” (NPO-30573), NASA Tech 
Briefs , Vol. 28, No. 10 (October 2004), 
page 36. Unlike those prior systems that 
depend on detection of sounds and/ or 
lights emitted by emergency vehicles, 
this system is not subject to severe range 
limitations. This system can be retrofit- 
ted into any pre-existing traffic-light- 
control system, without need to modify 
that system other than to make a mini- 
mal number of wire connections between 
the two systems. 

This system comprises several subsys- 
tems, including a transponder and inter- 
face circuitry on each emergency vehi- 
cle, a monitoring and control unit at 
each intersection equipped with traffic 
lights, and a wide-area two-way radio 
communication network that connects 
the emergency vehicles and intersection 
units. Computers in the various intersec- 
tions and vehicle units run special-pur- 


pose software that implements the traf- 
fic-light-preemption scheme. The opera- 
tions of the intersection and vehicle 
units are synchronized by use of Global 
Positioning System (GPS) timing signals. 
The transponder in each vehicle esti- 
mates its own position and velocity by 
use of GPS signals, deductive (“dead”) 
reckoning, data from the onboard diag- 
nostic (OBD) computer of the vehicle, 
and/or triangulation of beacon signals. 

When the operator of an emergency 
vehicle turns on its flashing lights and 
sirens in response to a request for an 
emergency response, the transponder 
unit goes into action, reading the OBD 
data to determine speed and accelera- 
tion, and reading and gathering further 
navigational data as described above. 
The position, velocity, and acceleration 
data are combined with vehicle-identifi- 
cation data in a prescribed format, and 
the resulting set of data is transmitted to 
the intersections within communication 
range of the transponder. 

In each intersection unit that receives 
such a data signal, a processor estimates 
the time of arrival of the vehicle, com- 
pares it with the estimated times of ar- 
rival of other emergency vehicles ap- 
proaching the intersection, and 
determines which vehicle will arrive first. 
The intersection unit notifies the trans- 
ponders of all emergency vehicles of a 


potential conflict and states, as part of 
the notice, which vehicle has the right of 
way. At the same time, the processor col- 
lects information on the current opera- 
tion of the traffic lights at the intersec- 
tion and calculates when pedestrians 
should be alerted not to cross and when 
preemption of the traffic lights should 
start. When preemption starts, the traffic 
lights are augmented by textual displays 
of a message that emergency vehicles are 
approaching and graphical displays indi- 
cating the direction (s) of approach. 
Once the emergency vehicles have 
passed through the intersection, normal 
operation of the traffic lights is resumed. 

This work was done by Conrad Foster and 
Aaron Bachelder of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP (see page 1 ). 

In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Innovative Technology Assets Management 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-30604, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 


Vehicle Transponder for Preemption of Traffic Lights 

This unit provides timely information on statuses of vehicles and intersections. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


The purpose of this article is to de- 
scribe, in more detail, the transponder 
installed in each vehicle that participates 
in the emergency traffic-light-preemp- 
tion system described in the immedi- 
ately preceding article. The transponder 
(see figure) is a fully autonomous data- 
collection, data-processing, information- 
display, and communication subsystem 
that performs robustly in preemption of 


traffic lights and monitoring of the sta- 
tuses of street intersections. 

This transponder monitors the condi- 
tion of the emergency vehicle in which it 
is installed and determines when the ve- 
hicle has been placed in an emergency- 
response condition with its siren and/ or 
warning lights activated. Upon detection 
of such a condition, the transponder col- 
lects real-time velocity and acceleration 


data from the onboard diagnostic 
(OBD) computer of the vehicle. For this 
purpose, the transponder contains an 
OBD interface circuit, including a mi- 
croprocessor that determines the manu- 
facturer and model of the vehicle and 
then sends the appropriate commands 
to the OBD computer requesting the 
speed and acceleration data. At the same 
time, data from an onboard navigation 
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An Emergency-Vehicle Transponder contains electronic circuits designed by NASA's Jet Propulsion Lab- 
oratory. The transponder is packaged such that it can be easily mounted in the vehicle in less than one 
hour. 


system are collected to determine the lo- 
cation and the heading of the vehicle. 
Then acceleration, speed, position, and 
heading data are processed and com- 
bined with a vehicle-identification num- 
ber and the resulting set of data is trans- 


mitted to monitoring and control units 
located at all intersections within com- 
munication range. 

When the unit at an intersection de- 
termines that this vehicle is approaching 
and has priority to preempt the intersec- 


tion, it transmits a signal declaring the 
priority and the preemption to all partic- 
ipating vehicles (including this one) in 
the vicinity. If the unit at the intersection 
has determined that other participating 
vehicles are also approaching the inter- 
section, then this unit also transmits, to 
the vehicle that has priority, a message 
that the other vehicles are approaching 
the same intersection. The texts of these 
messages, plus graphical symbols that 
show the directions and numbers of the 
approaching vehicles, are presented on 
the display panel of a computer that is 
part of the transponder. 

This work was done by Conrad Foster and 
Aaron Bachelder of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP (see page 1 ). 

In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Innovative Technology Assets Management 

jpl 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-30607, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 


Automated Announcements of Approaching Emergency Vehicles 

Pedestrians would be given advance warning. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Street intersections that are equipped 
with traffic lights would also be equipped 
with means for generating audible an- 
nouncements of approaching emergency 
vehicles, according to a proposal. The 
means to generate the announcements 
would be implemented in the intersec- 
tion-based subsystems of emergency traf- 
fic-light-preemption systems like those 
described in the two immediately preced- 
ing articles and in “Systems Would Pre- 
empt Traffic Lights for Emergency Vehi- 
cles” (NPO-30573), NASA Tech Briefs, Vol. 
28, No. 10 (October 2004), page 36. 

Preempting traffic lights is not, by it- 
self, sufficient to warn pedestrians at af- 
fected intersections that emergency ve- 
hicles are approaching. Automated 
visual displays that warn of approaching 
emergency vehicles can be helpful as a 
supplement to preemption of traffic 
lights, but experience teaches that for a 


variety of reasons, pedestrians often do 
not see such displays. Moreover, in noisy 
and crowded urban settings, the lights 
and sirens on emergency vehicles are 
often not noticed until a few seconds be- 
fore the vehicles arrive. 

According to the proposal, the traf- 
fic-light preemption subsystem at each 
intersection would generate an audible 
announcement — for example, “emer- 
gency vehicle approaching, please clear 
intersection” — whenever a preemp- 
tion was triggered. The subsystem 
would estimate the time of arrival of an 
approaching emergency vehicle by use 
of vehicle identity, position, and time 
data from one or more sources that 
could include units connected to traffic 
loops and/or transponders connected 
to diagnostic and navigation systems in 
participating emergency vehicles. The 
intersection-based subsystem would 


then start the announcement far 
enough in advance to enable pedestri- 
ans to leave the roadway before any 
emergency vehicles arrive. 

This work was done by Aaron Bachelder 
and Conrad Foster of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP (see page 1 ). 

In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Innovative Technology Assets Management 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-30623, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 
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Intersection Monitor for Traffic-Light-Preemption System 

This unit provides real-time phase data essential for effective preemption. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 



The Intersection Monitor, shown here with its covers off, provides real-time data on the phases of traf- 
fic lights at an intersection, without interfering with the traffic-light control circuitry. 


The figure shows an intersection mon- 
itor that is a key subsystem of an emer- 
gency traffic-light-preemption system 
that could be any of the systems de- 
scribed in the three immediately preced- 
ing articles and in “Systems Would Pre- 
empt Traffic Lights for Emergency 
Vehicles” (NPO-30573), NASA Tech Briefs, 
Vol. 28, No. 10 (October 2004), page 36. 
This unit is so named because it is in- 
stalled at an intersection, where it moni- 
tors the phases (in the sense of timing) 
of the traffic lights. The mode of opera- 
tion of this monitor is independent of 
the type of traffic-light-controller hard- 
ware or software in use at the intersec- 
tion. Moreover, the design of the moni- 
tor is such that (1) the monitor does not, 
by itself, affect the operation of the traf- 
fic-light controller and (2) in the event 
of a failure of the monitor, the traffic- 
light controller continues to function 
normally (albeit without preemption) . 

The monitor is installed in series with 
the traffic-light controller at an intersec- 
tion. The control signals of interest are 
monitored by use of high-impedance 
taps on affected control lines. These taps 
are fully isolated and further protected 
by high-voltage diodes that prevent any 
voltages or short circuits that arise within 
the monitor from affecting the con- 
troller. The signals from the taps are 
processed digitally and cleaned up by use 


of high-speed logic gates, and the result- 
ing data are passed on to other parts of 
the traffic-light-preemption intersection 
subsystem. The data are compared con- 
tinuously with data from vehicles and 
used to calculate timing for reliable pre- 
emption of the traffic lights. The pedes- 
trian crossing at the intersection is also 
monitored, and pedestrians are warned 
not to cross during preemption. 

This work was done by Aaron Bachelder 
and Conrad Foster of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP (see page 1 ). 


In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Innovative Technology Assets Management 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-30612, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 


Full-Duplex Digital Communication on a Single Laser Beam 

The laser beam would be transmitted witb one modulation and retroreflected with another modulation. 

Goddard Space Flight Center, Greenbelt, Maryland 


A proposed free-space optical communi- 
cation system would operate in a full-duplex 
mode, using a single constant-power laser 
beam for transmission and reception of bi- 
nary signals at both ends of the free-space 
optical path. The system was conceived for 
two-way data communication between a 
ground station and a spacecraft in a low 
orbit around the Earth. It has been esti- 
mated that in this application, a data rate of 
10 kb/s could be achieved at a ground-sta- 
tion-to-spacecraft distance of 320 km, using 
a laser power of only 100 mW. The basic 
system concept is also applicable to terres- 
trial free-space optical communications. 

The system (see figure) would include a 


diode laser at one end of the link (origi- 
nally, the ground station) and a liquid-crys- 
tal-based retroreflecting modulator at the 
other end of the link (originally, the space- 
craft). At the laser end, the beam to be 
transmitted would be made to pass through 
a quarter-wave plate, which would convert 
its linear polarization to right circular polar- 
ization. For transmission of data from the 
laser end to the retroreflector end, the laser 
beam would be modulated with subcarrier 
phase-shift keying (SC7PSK) . The transmit- 
ted beam would then pass through an aper- 
ture-sharing element (ASE) — basically, a 
mirror with a hole in it, used to separate the 
paths of the transmitted and received light 


beams. The transmitted beam would con- 
tinue outward through a telescope (which, 
in the original application, would be 
equipped with a spacecraft-tracking system) 
that would launch the transmitted beam 
along the free-space optical path to the 
retroreflector end. 

At the retroreflector end, a portion of 
the received laser beam would be sent to a 
demodulator for detection of the SC-PSK 
signal. For transmitting data to the laser 
end, the retroreflected portion of the re- 
ceived laser beam would be modulated 
with circular-polarization keying (CPK), in 
which left circular polarization signifies a bi- 
nary level (“1” in this case) and right circu- 
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The Laser at One End of the free-space optical path would provide all of the beam power needed for transmission of data signals in both directions along the path. 


lar polarization signifies the other binary 
level (“0” in this case). Hence, to transmit 
“0,” the retroreflecting modulator would 
leave the right circular polarization of the 
retroreflected beam unchanged; to trans- 
mit “1,” the retroreflecting modulator 
would flip the polarization of the reflected 
beam to left circular. Full-duplex operation 
would be possible because the CPK and the 
SC-PSK would be transparent to each other. 

At the laser end, the reflected, CPK- 
modulated beam would return through 
the telescope and would then be reflected 
by the ASE into a receiver subsystem. A 
beam splitter would divert 0.2 percent of 
the beam power to a camera in the track- 
ing system. The remainder of the beam 
would pass through the beam splitter to a 
quarter-wave plate, which would convert 
the circular polarization to two orthogonal 
linear polarizations. A polarizing beam 
splitter would then split the light in these 
two polarizations so that photons corre- 


sponding to “0” would go to one photode- 
tector and photons corresponding to “1” 
would go to another photodetector. 

It should be emphasized that this arrange- 
ment would yield a nonzero photodetector 
output of nominally the same magnitude for 
either “0” or “1.” This is fundamentally dif- 
ferent from on-off keying (OOK) , in which 
“0” or “1” is represented by the absence or 
presence, respectively, of a signal. Taking 
advantage of this, prior to final digitization 
of the return signal at “0” or “1,” the out- 
put of the “0” photodetector could be in- 
verted, then subtracted from the output of 
the “1” photodetector to obtain twice the 
signal-to-noise ratio achievable in OOK 

The receiver subsystem would include 
Faraday-anomalous-dispersion optical fil- 
ters (FADOFs), which would reject back- 
ground light to such a high degree that 
the system could operate over a long path 
during daytime. The FADOFs would es- 
sentially prevent skylight from reaching 


the photodetectors while allowing about 
80 percent of the signal photons to pass 
through. Without the FADOFs, it would 
be necessary to increase the laser power 
by a factor of 10 for daytime operation. 

This work was done by D. A. Hazzard, J. A. 
MacCannell, G. Lee, E. R Selves, D. Moore, J. 
A. Payne, C. D. Garrett, N. Dahlstrom, and T. 
M. Shay of New Mexico State University for 
Goddard Space Flight Center. Further infor- 
mation is contained in a TSP (seepage 1 ). 

In accordance with Public Law 96-51 7, the 
contractor has elected to retain title to this inven- 
tion. Inquiries concerning rights for its commer- 
cial use should be addressed to: 

New Mexico State University Las Cruces 
P.O.Box 30001 
Las Cruces, NM 88003 
Refer to GSC-14759-1, volume and number 
of this NASA Tech Briefs issue, and the page 
number. 


Stabilizing Microwave Frequency of a Photonic Oscillator 

Microwave frequency is stabilized by stabilizing optical frequency to an atomic transition. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A scheme for stabilizing the frequency 
of a microwave signal is proposed that ex- 
ploits the operational characteristics of a 
coupled optoelectronic oscillator (COEO) 
and related optoelectronic equipment. An 
essential element in the scheme is a fiber 
mode-locked laser (MLL), the optical fre- 
quency of which is locked to an atomic 
transition. In this scheme, the optical fre- 
quency stability of the mode-locked laser is 
transferred to that of the microwave in the 


same device. Relative to prior schemes for 
using wideband optical frequency comb to 
stabilize microwave signals, this scheme is 
simpler and lends itself more readily to im- 
plementation in relatively compact, 
rugged equipment. The anticipated devel- 
opment of small, low-power, lightweight, 
highly stable microwave oscillators based 
on this scheme would afford great benefits 
in communication, navigation, metrology, 
and fundamental sciences. 


COEOs of various designs, at various 
stages of development, in some cases 
called by different names, have been de- 
scribed in a number of prior NASA Tech 
Briefs articles. A COEO is an optoelec- 
tronic apparatus that generates both 
short (picosecond) optical pulses and a 
steady microwave signal having an ultra- 
high degree of spectral purity. The term 
“coupled optoelectronic” in the full 
name of such an apparatus signifies that 
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The COEO Generates a Comblike Optical Spectrum that is used to generate a microwave signal. The 
frequency stability of a diode laser signal locked to the frequency of an atomic transition is transferred 
to the microwave signal. 


its optical and electronic oscillations are 
coupled to each other in a single device. 

The present frequency-stabilization 
scheme is best described indirecdy by de- 
scribing the laboratory apparatus used to 
demonstrate it. The apparatus (see figure) 
includes a COEO that generates a comb- 
like optical spectrum, the various fre- 
quency components of which interfere, 
producing short optical pulses. This spec- 
trum is centered at a nominal wavelength 
of 1,560 nm. The spectrum separation of 
this comb is about 10 GHz, as determined 
primarily by the length of an optical loop 
and the bandpass filter in the microwave 
feedback loop. The optical loop serves as 
microwave resonator having a very high 
value of the resonance quality factor ( 0 . 
The optical frequency of MLL is then sta- 
bilized by locking it to an atomic transition 
as described below. 

The COEO contains a tunable 1-nm 
band-pass optical filter and a piezoelectric- 
transducer (PZT) drum over which a 
stretch of fiber is wound. The 1-nm-wide 
pass band of the filter provides coarse tun- 
ing to overlap the frequency comb with the 
atomic transition frequency. Controlled 
stretching of the fiber by means of the PZT 
drum can be used in conjunction with tem- 
perature control for locking the laser fre- 
quency. To reference to an atomic reso- 
nance at 780 nm in this demonstration 
setup, the optical output of the COEO at 
1,560 nm is fed through an erbium-doped- 
fiber amplifier (EDFA) to a frequency dou- 
bler in the form of a periodically poled 


lithium niobate (PPLN) crystal. The fre- 
quency-doubled output is combined with 
the output of a separate frequency-stabi- 
lized diode laser at a photodetector. As de- 
scribed thus far, the two 780-nm laser sub- 
systems are nominally independent of each 
other and can, therefore, operate at differ- 
ent frequencies. Hence, at the photodetec- 
tor, the two laser beams interfere, so that 
the output of the photodetector includes a 
beat note (a component at the difference 
between the two laser frequencies) . 

The beat note is used to stabilize the rel- 
ative frequency between the two optical 
signals through a simple electronic feed- 
back subsystem that adjusts the voltage ap- 
plied to the PZT to lock the optical fre- 
quency of the COEO to that of the diode 
laser. The diode laser is frequency stabi- 


lized to the atomic absorption of Rb vapor 
through frequency modulation (FM) satu- 
ration spectroscopy. The fractional fre- 
quency stability of it has been shown to be 
10 -12 at 1 second. After further optimiza- 
tion of design to minimize destabilizing ef- 
fects, it may be possible to attain a long- 
term stability at 10 -13 . Such optical 
frequency stability can be transferred to 
the microwave in the COEO device where 
the optical and microwave oscillators are 
coupled, and hence producing a highly 
stable microwave signal. 

This work was done by Lute Maleki, Nan 
Yu, and Meirong Tu of Caltech for NASA’s J et 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP (seepage 1 ). 
NPO-43026 


Microwave Oscillators Based on Nonlinear WGM Resonators 

Optical signals are phase-modulated with spectrally pure microwave signals. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Optical oscillators that exploit reso- 
nantly enhanced four-wave mixing in non- 
linear whispering-gallery-mode (WGM) 
resonators are under investigation for po- 
tential utility as low-power, ultra-miniature 
sources of stable, spectrally pure mi- 
crowave signals. There are numerous po- 
tential uses for such oscillators in radar 
systems, communication systems, and sci- 
entific instrumentation. 

The resonator in an oscillator of this 
type is made of a crystalline material that 
exhibits cubic Kerr nonlinearity, which 
supports the four-photon parametric 
process also known as four-wave mixing. 
The oscillator can be characterized as all- 
optical in the sense that the entire process 
of generation of the microwave signal 


takes place within the WGM resonator. 
The resonandy enhanced four-wave mix- 
ing yields coherent, phase-modulated op- 
tical signals at frequencies governed by the 
resonator structure. The frequency of the 
phase-modulation signal, which is in the 
microwave range, equals the difference 
between the frequencies of the optical sig- 
nals; hence, this frequency is also gov- 
erned by the resonator structure. Hence, 
further, the microwave signal is stable and 
can be used as a reference signal. 

The figure schematically depicts the 
apparatus used in a proof-of-principle ex- 
periment. Linearly polarized pump light 
was generated by an yttrium aluminum 
garnet laser at a wavelength of 1.32 pm. 
By use of a 90:10 fiber-optic splitter and 


optical fibers, some of the laser light was 
sent into a delay line and some was trans- 
mitted to one face of glass coupling 
prism, that, in turn, coupled the laser 
light into a crystalline CaF 2 WGM disk 
resonator that had a resonance quality 
factor (Q) of 6 x 10 9 . The output light of 
the resonator was collected via another 
face of the coupling prism and a single- 
mode optical fiber, which transmitted 
the light to a 50:50 fiber-optic splitter. 
One output of this splitter was sent to a 
slow photodiode to obtain a DC signal 
for locking the laser to a particular res- 
onator mode. The other output of this 
splitter was combined with the delayed 
laser signal in another 50:50 fiber-optic 
splitter used as a combiner. The output 
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This Laboratory Setup was used to demonstrate the basic principle of generating a microwave signal 
in a low-threshold hyperparametric process. 


of the combiner was fed to a fast photo- 
diode that demodulated light and gener- 
ated microwave signal. 

In this optical configuration, the res- 
onator was incorporated into one arm of 
a Mach-Zehnder interferometer, which 
was necessary for the following reasons: 
It was found that when the output of the 
resonator was sent directly to a fast pho- 
todiode, the output of the photodiode 


did not include a measurable microwave 
signal. However, when the resonator was 
placed in an arm of the interferometer 
and the delay in the other arm was set at 
the correct value, the microwave signal 
appeared. Such behavior is distinctly 
characteristic of phase-modulated light. 

The phase-modulation signal had a fre- 
quency of about 8 GHz, corresponding to 
the free spectral range of the resonator. 


The spectral width of this microwave signal 
was less than 200 Hz. The threshold pump 
power for generating the microwave signal 
was about 1 mW. It would be possible to re- 
duce the threshold power by several orders 
of magnitude if resonators could be made 
from crystalline materials in dimensions 
comparable to those of microresonators 
heretofore made from fused silica. 

This work was done by Lute Maleki, Au- 
drey Matsko, Anatoliy Savchenkov, and 
Dmitry Strekalov of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP ( see page 1 ). 

In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Innovative Technology Assets Management 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-41074, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 


Pointing Reference Scheme for Free-Space Optical 
Communications Systems 

A technique is proposed for referencing infrared transmit lasers with silicon detectors. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A scheme is proposed for referencing 
the propagation direction of the transmit 
laser signal in pointing a free-space optical 
communications terminal. This recendy 
developed scheme enables the use of low- 
cost, commercial silicon-based sensors for 
tracking the direction of the transmit laser, 
regardless of the transmit wavelength. 
Compared with previous methods, the 
scheme offers some advantages of less me- 
chanical and optical complexity and avoids 


expensive and exotic sensor technologies. 

In free-space optical communications, 
the transmit beam must be accurately 
pointed toward the receiver in order to 
maintain the communication link. The 
current approaches to achieve this func- 
tion call for part of the transmit beam to 
be split off and projected onto an optical 
sensor used to infer the pointed direction. 
This requires that the optical sensor be 
sensitive to the wavelength of the transmit 


laser. If a different transmit wavelength is 
desired, for example to obtain a source 
capable of higher data rates, this can be- 
come quite impractical because of the un- 
availability or inefficiency of sensors at 
these wavelengths. The innovation pro- 
posed here decouples this requirement by 
allowing any transmit wavelength to be 
used with any sensor. 

We have applied this idea to a particular 
system that transmits at the standard 



An Example Implementation is shown of the much simpler pointing reference scheme. 
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telecommunication wavelength of 1,550 
nm and uses a silicon-based sensor, sensi- 
tive from 0.5 to 1.0 micrometers, to deter- 
mine the pointing direction. The scheme 
shown in the figure involves integrating a 
low-power 980-nm reference or boresight 
laser beam coupled to the 1,550-nm trans- 
mit beam via a wavelength-division-multi- 
plexed fiber coupler. Both of these signals 
propagate through the optical fiber where 
they achieve an extremely high level of co- 
alignment before they are launched into 
the telescope. The telescope uses a 
dichroic beam splitter to reflect the 980- 
nm beam onto the silicon image sensor (a 
quad detector, charge-coupled device, or 
active-pixel-sensor array) while the 1,550- 
nm signal beam is transmitted through the 
optical assembly toward the remotely lo- 
cated receiver. Since the 980-nm reference 
signal originates from the same single- 
mode fiber-coupled source as the transmit 
signal, its position on the sensor is used to 


accurately determine the propagation di- 
rection of the transmit signal. 

The optics are considerably simpler in 
the proposed scheme due to the use of a 
single aperture for transmitting and re- 
ceiving. Moreover, the issue of mechani- 
cal misalignment does not arise because 
the reference signal and transmitted 
laser beams are inherently co-aligned. 
The beam quality of the 980-nm refer- 
ence signal used for tracking is required 
to be circularly symmetric and stable at 
the tracking-plane sensor array in order 
to minimize error in the centroiding al- 
gorithm of the pointing system. How- 
ever, since the transmit signal is delivered 
through a fiber that supports a single 
mode at 1,550 nm, propagation of 
higher order 980-nm modes is possible. 
Preliminary analysis shows that the over- 
all mode profile is dominated by the fun- 
damental mode, giving a near symmetric 
profile. The instability of the mode was 


also measured and found to be negligi- 
ble in comparison to the other error con- 
tributions in the centroid position on the 
sensor array. 

This work was done by Malcolm Wright, 
Gerardo Ortiz, and Muthu Jeganathan of 
Caltech for NASA’s Jet Propulsion Labo- 
ratory. Further information is contained in 
a TSP (see page 1 ). 

In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Innovative Technology Assets Management 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-30606, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 


High-Level Performance Modeling of SAR Systems 

NASA’s Jet Propulsion Laboratory , Pasadena , California 


SAUSAGE (Still Another Utility for 
SAR Analysis that’s General and Extensi- 
ble) is a computer program for modeling 
(see figure) the performance of syn- 


thetic-aperture radar (SAR) or interfero- 
metric synthetic-aperture radar (InSAR 
or 1FSAR) systems. The user is assumed 
to be familiar with the basic principles of 


SAR imaging and interferometry. Given 
design parameters (e.g., altitude, power, 
and bandwidth) that characterize a radar 
system, the software predicts various per- 
formance metrics (e.g., signal-to-noise 
ratio and resolution). SAUSAGE is in- 
tended to be a general software tool for 
quick, high-level evaluation of radar de- 
signs; it is not meant to capture all the 
subtleties, nuances, and particulars of 
specific systems. SAUSAGE was written to 
facilitate the exploration of engineering 
tradeoffs within the multidimensional 
space of design parameters. Typically, 
this space is examined through an itera- 
tive process of adjusting the values of the 
design parameters and examining the ef- 
fects of the adjustments on the overall 
performance of the system at each itera- 
tion. The software is designed to be mod- 
ular and extensible to enable considera- 
tion of a variety of operating modes and 
antenna beam patterns, including, for 
example, strip-map and spotlight SAR ac- 
quisitions, polarimetry, burst modes, and 
squinted geometries. 

This program was written by Curtis Chen 
of Caltech for NASA’s Jet Propulsion Lab- 
oratory. Further information is contained in 
a TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-43373. 


System 



Logical Block Diagram shows how a radar system is modeled by SAUSAGE. Blocks represent software 
objects, which describe physical aspects of the system. Arrays are indicated by stacked blocks. Arrows 
represent pointers to other objects. 
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Spectral Analysis Tool 6.2 for Windows 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Spectral Analysis Tool 6.2 is the latest 
version of a computer program that as- 
sists in analysis of interference between 
radio signals of the types most commonly 
used in Earth/spacecraft radio commu- 
nications. [An earlier version was re- 
ported in “Software for Analyzing 
Earth/ Spacecraft Radio Interference” 
(NPO-20422), NASA Tech Briefs, Vol. 25, 
No. 4 (April 2001), page 52.] SAT 6.2 cal- 
culates signal spectra, bandwidths, and 
interference effects for several families of 
modulation schemes. Several types of fil- 
ters can be modeled, and the program 


calculates and displays signal spectra 
after filtering by any of the modeled fil- 
ters. The program accommodates two si- 
multaneous signals: a desired signal and 
an interferer. The interference-to-signal 
power ratio can be calculated for the fil- 
tered desired and interfering signals. 
Bandwidth-occupancy and link-budget 
calculators are included for the user’s 
convenience. SAT 6.2 has a new software 
structure and provides a new user inter- 
face that is both intuitive and conven- 
ient. SAT 6.2 incorporates multi-tasking, 
multithreaded execution, virtual mem- 


ory management, and a dynamic link li- 
brary. SAT 6.2 is designed for use on 32- 
bit computers employing Microsoft Win- 
dows operating systems. 

This program was written by Feiming Mor- 
gan, Miles Sue, Ted Peng, Harry Tan, and 
Robert Liang of Caltech and Peter Kinman of 
California State University, Fresno, for NASA’s 
Jet Propulsion Laboratory. Further informa- 
tion is contained in a TSP ( seepage 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-43129. 


Multi-Platform Avionics Simulator 

NASA’s Jet Propulsion Laboratory, Pasadena, California 



Multi-Platform Avionics Simulator 
(MPAvSim) is a software library for de- 
velopment of simulations of avionic 
hardware. MPAvSim facilitates simula- 
tion of interactions between flight soft- 
ware and such avionic peripheral 
equipment as telecommunication de- 
vices, thrusters, pyrotechnic devices, 
motor controllers, and scientific instru- 
ments. MPAvSim focuses on the behav- 
ior of avionics as seen by flight soft- 
ware, rather than on performing 
high-fidelity simulations of dynamics. 
However, MPAvSim is easily integrable 
with other programs that do perform 
such simulations. MPAvSim makes it 
possible to do real-time partial hard- 
ware-in-the-loop simulations. An 
MPAvSim simulation consists of execu- 
tion chains (see figure) represented by 
flow graphs of models, defined here as 
stateless procedures that do some 
work. During a simulation, MPAvSim 
walks the execution chain, running 
each model in turn. Using MPAvSim, 
flight software can be run against a 
spacecraft that is all simulation, all 
hardware, or part hardware and part 


simulation. With respect to a specific 
piece of hardware, either the hardware 
itself or its simulation can be plugged 
in without affecting the rest of the sys- 
tem. Thus, flight software can be tested 
before hardware is available, and as 
items of hardware become available, 
they can be substituted for their simu- 
lations, with minimal disruption. 


This program was written by Micah Clark, 
Robert Steinke, and Elihu McMahon of Cal- 
tech for NASA’s Jet Propulsion Labora- 
tory. Further information is contained in a 
TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41277. 
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Silicon-Based Optical Modulator With Ferroelectric Layer 

This device would remain switched even with power turned off. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


According to a proposal, a silicon diox- 
ide layer in a high-speed, low-power, sili- 
con-based electro-optical modulator 
would be replaced by a layer of lead zir- 
conate titanate or other ferroelectric 
oxide material. The purpose of this mod- 
ification is to enhance the power perform- 
ance and functionality of the modulator. 

In its unmodified form, the particular sil- 
icon-based electro-optical modulator is of 
an advanced design that overcomes the 
speed limitation of prior silicon-based elec- 
tro-optical modulators. Whereas modula- 
tion frequencies of such devices had been 
limited to about 20 MHz, this modulator 
can operate at modulation frequencies as 
high as 1 GHz. This modulator can be char- 
acterized as a silicon-waveguide-based 


metal oxide/ semiconductor (MOS) capac- 
itor phase shifter in which modulation of 
the index of refraction in silicon is obtained 
by exploiting the fr ee-charge-carrier-plasma 
dispersion effect. As shown in the figure, 
the modulator includes an n-doped crys- 
talline silicon slab (the silicon layer of a sili- 
con-on-insulator wafer) and a p-doped 
polycrystalline silicon rib with a gate oxide 
layer (the aforementioned silicon dioxide 
layer) sandwiched between them. 

Under accumulation conditions, the 
majority charge carriers in the silicon 
waveguide modify the index of refraction 
so that a phase shift is induced in the op- 
tical mode propagating in the waveguide. 
The advantage of using an MOS capacitor 
phase shifter is that it is possible to 


achieve high modulation speed because 
there are no slow carrier-generation 
or -recombination processes involved in 
the accumulation operation. 

The main advantage of the proposed 
substitution of a ferroelectric oxide layer 
for the silicon dioxide layer would arise 
from the spontaneous polarization effect 
of the ferroelectric layer: This sponta- 
neous polarization would maintain accu- 
mulation conditions in the absence of 
applied voltage. Consequently, once the 
device had been switched to a given opti- 
cal state, it would remain in that state, 
even in the absence of applied voltage 
(in other words, even with power turned 
off). A secondary advantage is that be- 
cause the ferroelectric layer would have 
an index of refraction larger than that of 
silicon dioxide, there could be some re- 
duction of optical losses attributable to 
fabrication of the modulator. 

This work was done by Douglas Sheldon of 
Caltech for NASA’s Jet Propulsion Labo- 
ratory. Further information is contained in 
a TSP (see page 1 ). 

In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Innovative Technology Assets Management 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-40935, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 



A Ferroelectric Oxide Layer would be substituted for the gate oxide layer in this silicon-based electro- 
optical modulator. 


Multiplexing Transducers Based on Tunnel-Diode Oscillators 

Compact, low-power transducers could operate over wide temperature ranges. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Multiplexing and differential transduc- 
ers based on tunnel-diode oscillators 
(TDOs) would be developed, according to 
a proposal, for operation at very low 
and/or widely varying temperatures in ap- 
plications that involve requirements to 
minimize the power and mass of trans- 
ducer electronic circuitry. It has been 
known since 1975 that TDOs are useful for 


making high-resolution (of the order of 
10 -9 ) measurements at low temperatures. 
Since that time, TDO transducers have 
been found to offer the following addi- 
tional advantages, which the present pro- 
posal is intended to exploit: 

• TDO transducers can operate at temper- 
atures ranging from 1 K to about 400 K 
Most electronic components other than 


tunnel diodes do not operate over such 
a wide temperature range. 

• TDO transducers can be made to oper- 
ate at very low power — typically, <1 mW. 

• Inasmuch as the response of a TDO 
transducer is a small change in an arbi- 
trarily set oscillation frequency, the 
outputs of many TDOs operating at 
sufficiently different set frequencies 
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can be multiplexed through a single 
wire. 

• Inasmuch as frequencies can be easily 
subtracted by means of mixing circuitry, 
one can easily use two TDOs to make dif- 
ferential measurements. Differential 
measurements are generally more pre- 
cise and less susceptible to environ- 
mental variations than are absolute 
measurements. 

• TDO transducers are tolerant to ioniz- 
ing radiation. 

• Ultimately, the response of a TDO 
transducer is measured by use of a fre- 
quency counter. Because frequency 
counting can be easily implemented by 
use of clock signals available from most 


microprocessors, it is not necessary to 
incorporate additional readout cir- 
cuitry that would, if included, add to 
the mass and power consumption of 
the transducer circuitry. 

In one example of many potential vari- 
ations on the basic theme of the pro- 
posal, the figure schematically depicts a 
conceptual differential-pressure trans- 
ducer containing a symmetrical pair of 
TDOs. The differential pressure would be 
exerted on an electrically conductive and 
grounded diaphragm, which, at zero dif- 
ferential pressure, would nominally be 
sprung to a middle position between two 
capacitor plates that would be parts of the 
two TDOs. The frequencies of the two 


TDOs would vary in opposite directions 
as variations in differential pressure bent 
the diaphragm away from one capacitor 
plate and toward the other. The outputs 
of the TDOs would be mixed and low- 
pass filtered to obtain a signal at the dif- 
ference between the frequencies of the 
two TDOs. The difference frequency 
would be measured by a frequency 
counter and converted to differential 
pressure by a computer. 

This work was done by Talso Chui, Kon- 
stantin Penanen, and Joseph Young of Cal- 
tech for NASA’s Jet Propulsion Labora- 
tory. Further information is contained in a 
TSP (see page 1 ). 

NPO-43079 
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Q Scheduling With Automatic 
Resolution of Conflicts 

DSN Requirement Scheduler is a com- 
puter program that automatically sched- 
ules, reschedules, and resolves conflicts 
for allocations of resources of NASA’s 
Deep Space Network (DSN) on the basis 
of ever-changing project requirements 
for DSN services. As used here, “re- 
sources” signifies, primarily, DSN anten- 
nas, ancillary equipment, and times dur- 
ing which they are available. Examples 
of project-required DSN services include 
arraying, segmentation, very-long-base- 
line interferometry, and multiple space- 
craft per aperture. Requirements can in- 
clude periodic reservations of specific or 
optional resources during specific time 
intervals or within ranges specified in 
terms of starting times and durations. 
This program is built on the Automated 
Scheduling and Planning Environment 
(ASPEN) software system (aspects of 
which have been described in previous 
NASA Tech Briefs articles), with cus- 
tomization to reflect requirements and 
constraints involved in allocation of 
DSN resources. Unlike prior DSN-re- 
source-scheduling programs that make 
single passes through the requirements 
and require human intervention to re- 
solve conflicts, this program makes re- 
peated passes in a continuing search for 
all possible allocations, provides a best- 
effort solution at any time, and presents 
alternative solutions among which users 
can choose. 

This program was written by Bradley 
Clement and Steve Schaffer of Caltech for 

NASA’s Jet Propulsion Laboratory. Fur- 
ther information is contained in a TSP (see 
page 1). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41313. 


Symbolic Constraint 
Maintenance Grid 

Version 3.1 of Symbolic Constraint 
Maintenance Grid (SCMG) is a software 
system that provides a general concep- 
tual framework for utilizing pre-existing 
programming techniques to perform 
symbolic transformations of data. SCMG 
also provides a language (and an associ- 
ated communication method and proto- 


col) for representing constraints on the 
original non-symbolic data. SCMG pro- 
vides a facility for exchanging informa- 
tion between numeric and symbolic com- 
ponents without knowing the details of 
the components themselves. In essence, 
it integrates symbolic software tools (for 
diagnosis, prognosis, and planning) with 
non-artificial-intelligence software . 
SCMG executes a process of symbolic 
summarization and monitoring of con- 
tinuous time series data that are being 
abstractly represented as symbolic tem- 
plates of information exchange. This 
summarization process enables such sym- 
bolic-reasoning computing systems as ar- 
tificial-intelligence planning systems to 
evaluate the significance and effects of 
channels of data more efficiently than 
would otherwise be possible. As a result 
of the increased efficiency in representa- 
tion, reasoning software can monitor 
more channels and is thus able to per- 
form monitoring and control functions 
more effectively. 

This work was done by Mark James of Cal- 
tech for NASA’s Jet Propulsion Labora- 
tory. Further information is contained in a 
TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-42001. 


3 Discerning Trends in Per- 
formance Across Multiple 
Events 

Mass Data is a computer program that 
enables rapid, easy discernment of 
trends in performance data across multi- 
ple flights and ground tests. The pro- 
gram can perform Fourier analysis and 
other functions for the purposes of fre- 
quency analysis and trending of all vari- 
ables. These functions facilitate identifi- 
cation of past use of diagnosed systems 
and of anomalies in such systems, and 
enable rapid assessment of related cur- 
rent problems. Many variables, for com- 
putation of which it is usually necessary 
to perform extensive manual manipula- 
tion of raw downlist data, are automati- 
cally computed and made available to all 
users, regularly eliminating the need for 
what would otherwise be an extensive 
amount of engineering analysis. Data 
from flight, ground test, and simulation 
are preprocessed and stored in one cen- 


tral location for instantaneous access 
and comparison for diagnostic and 
trending purposes. Rules are created so 
that an event log is created for every 
flight, making it easy to locate informa- 
tion on similar maneuvers across many 
flights. The same rules can be created 
for test sets and simulations, and are 
searchable, so that information on like 
events is easily accessible. 

This program was written by Simon Slater, 
Mike Hiltz, and Craig Rice of MacDonald 
Dettwiler Space Robotics for Johnson Space 
Center. For further information, contact the 
Johnson Technology Transfer Office at (281) 
483-3809. 

MSC-23824 


Magnetic Field Solver 

The Magnetic Field Solver computer 
program calculates the magnetic field 
generated by a group of collinear, cylin- 
drical axisymmetric electromagnet coils. 
Given the current flowing in, and the 
number of turns, axial position, and axial 
and radial dimensions of each coil, the 
program calculates matrix coefficients 
for a finite-difference system of equations 
that approximates a two-dimensional par- 
tial differential equation for the magnetic 
potential contributed by the coil. The 
program iteratively solves these finite-dif- 
ference equations by use of the modified 
incomplete Cholesky preconditioned- 
conjugate-gradient method. The total 
magnetic potential as a function of axial 
(z) and radial (r) position is then calcu- 
lated as a sum of the magnetic potentials 
of the individual coils, using a high-accu- 
racy interpolation scheme. Then the r 
and z components of the magnetic field 
as functions of rand z are calculated from 
the total magnetic potential by use of a 
high-accuracy finite-difference scheme. 
Notably, for the finite-difference calcula- 
tions, the program generates nonuni- 
form two-dimensional computational 
meshes from nonuniform one-dimen- 
sional meshes. Each mesh is generated in 
such a way as to minimize the numerical 
error for a benchmark one-dimensional 
magnetostatic problem. 

This program was written by Andrew V. 
Ilin of Muniz Engineering, Inc. for Johnson 
Space Center. For further information, con- 
tact the Johnson Technology Transfer Office 
at (281) 483-3809. 

MSC-23618 
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J Computing for Aiming a 
Spaceborne Bistatic-Radar 
Transmitter 

BISTAT is a computer program for 
use in aiming a spaceborne bistatic- 
radar transmitting antenna at a remote 
planet that has an atmosphere, such that 
after refraction by the atmosphere and 
reflection from the surface of the 
planet, the radar signal travels toward a 
receiver on Earth. BISTAT includes an 


algorithm that neglects atmospheric re- 
fraction and calculates a specular-reflec- 
tion point for a spacecraft at a given lo- 
cation. The specular-reflection point is 
then used as an initial guess for a modi- 
fied limb-track algorithm that takes at- 
mospheric refraction into account. The 
output of BISTAT for all spacecraft posi- 
tions of interest constitutes a pointing 
profile; the output data are in the form 
of an inertial-vector file and a Doppler- 
residual file. The inertial-vector file is 


used to command the attitude of the 
spacecraft; the Doppler-residual file is 
used to determine a downlink frequency 
file for the receiver. 

This program was written by Nicole Rap- 
paport of Caltech for NASA’s Jet Propul- 
sion Laboratory. Further information is 
contained in a TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41518. 
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$ 4-Vinyl-l ,.'!-I)ioxolane-2-()iie as an Additive for Li-Ion Cells 

Low-temperature charge/discharge capacity is increased. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 



The Proportion of Room-Temperature Capacity retained at -20 °C was determined, for each of two 
cells, from data obtained in the first five cycles of charge/discharge testing. 


Electrolyte additive 4-vinyl-l,3-diox- 
olane-2-one has been found to be prom- 
ising for rechargeable lithium-ion elec- 
trochemical cells. This and other 
additives, along with advanced elec- 
trolytes comprising solutions of LiPF 6 in 
various mixtures of carbonate solvents, 
have been investigated in a continuing 
effort to improve the performances of 
rechargeable lithium-ion electrochemi- 
cal cells, especially at low temperatures. 
In contrast to work by other researchers 
who have investigated the use of this ad- 
ditive to improve the high-temperature 
resilience of Li-ion cells, the current 
work involves the incorporation of 4- 
vinyl-l,3-dioxolane-2-one into quater- 
nary carbonate electrolyte mixtures, pre- 
viously optimized for low-temperature 
applications, resulting in improved low- 
temperature performance. 

The benefit afforded by 4-vinyl-l,3- 
dioxolane-2-one can be better under- 
stood in the light of relevant informa- 
tion from a number of prior NASA Tech 
Briefs articles about electrolytes and ad- 
ditives for such cells. To recapitulate: 
The loss of performance with decreasing 
temperature is attributable largely to a 
decrease of ionic conductivity and the 
increase in viscosity of the electrolyte. 
What is needed to extend the lower limit 
of operating temperature is a stable elec- 
trolyte solution with relatively small low- 
temperature viscosity, a large electric 
permittivity, adequate coordination be- 
havior, and appropriate ranges of solu- 
bilities of liquid and salt constituents. 
Whether the anode is made of graphitic 
or non-graphitic carbon, a film on the 
surface of the anode acts as a solid/elec- 
trolyte interface (SEt), the nature of 
which is critical to low-temperature per- 
formance. Desirably, the surface film 
should exert a chemically protective 
(passivating) effect on both the anode 
and the electrolyte, yet should remain 
conductive to lithium ions to facilitate 
intercalation and deintercalation of the 
ions into and out of the carbon during 
discharging and charging, respectively. 

The additives investigated previously 
include alkyl pyrocarbonates. Those ad- 


ditives help to improve low-temperature 
performances by giving rise to the for- 
mation of SEIs having desired proper- 
ties. The formation of the SEIs is be- 
lieved to be facilitated by products (e.g., 
C0 2 ) of the decomposition of these ad- 
ditives. These decomposition products 
are believed to react to form Li 2 C0 3 - 
based films on the carbon electrodes. 

The present additive, 4-vinyl-l,3-diox- 
olane-2-one, also helps to improve low- 
temperature performance by contribut- 
ing to the formation of SEIs having 
desired properties, but probably in a dif- 
ferent manner: It is believed that, as part 
of the decomposition process, the com- 
pound polymerizes on the surfaces of 
carbon electrodes. 

The effectiveness of 4-vinyl- 1,3-diox- 
olane-2-one as a performance-improving 
additive in several different LiPF 6 /car- 
bonate-solvent-mixture electrolytes was 
investigated in a series of charge/dis- 
charge tests of rechargeable lithium-ion 
electrochemical cells containing the 
electrolytes, at room temperature and at 


a temperature of -20 °C. In one pair of 
tests, the electrolyte comprised LiPF 6 
dissolved at a concentration of 1.0 M in 
a solvent mixture comprising 1 volume 
part of ethylene carbonate + 1 volume 
part of diethyl carbonate + 1 volume 
part of dimethyl carbonate + 2 volume 
parts of ethyl methyl carbonate. In one 
of the tests, no additive was included; in 
the other test, the electrolyte included 4- 
vinyl-l,3-dioxolane-2-one in the propor- 
tion of 1.5 weight percent. The results of 
these tests showed that the additive en- 
hanced low-temperature performance: 
specifically, the results showed that at 
-20 °C and at any given voltage, the cell 
containing the additive retained a 
greater proportion of its room-tempera- 
ture capacity than did the cell that did 
not contain the additive (see figure) . 

This work was done by Marshall Smart 
and Ratnakumar Bugga of Caltech for 
NASA’s Jet Propulsion Laboratory. Fur- 
ther information is contained in a TSP ( see 
page 1). 

NPO-40969 
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$ Probabilistic Prediction of Lifetimes of Ceramic Parts 

John H. Glenn Research Center, Cleveland, Ohio 


ANSYS/CARES/PDS is a software sys- 
tem that combines the ANSYS Probabilis- 
tic Design System (PDS) software with a 
modified version of the Ceramics Analysis 
and Reliability Evaluation of Structures 
Life (CARES/Life) Version 6.0 software. 
[A prior version of CARES/Life was re- 
ported in “Program for Evaluation of Reli- 
ability of Ceramic Parts” (LEW-16018), 
NASA Tech Briefs, Vol. 20, No. 3 (March 
1996), page 28.] CARES/Life models ef- 
fects of stochastic strength, slow crack 
growth, and stress distribution on the over- 
all reliability of a ceramic component. The 
essence of the enhancement in 
CARES/Life 6.0 is the capability to predict 


the probability of failure using results 
from transient finite-element analysis. 
ANSYS PDS models the effects of uncer- 
tainty in material properties, dimensions, 
and loading on the stress distribution and 
deformation. ANSYS/CARES/PDS ac- 
counts for the effects of probabilistic 
strength, probabilistic loads, probabilistic 
material properties, and probabilistic tol- 
erances on the lifetime and reliability of 
the component. Even failure probability 
becomes a stochastic quantity that can be 
tracked as a response variable. 
ANSYS/CARES/PDS enables tracking of 
all stochastic quantities in the design 
space, thereby enabling more precise 


probabilistic prediction of lifetimes of ce- 
ramic components. 

This work was done by Noel N. Nemeth and 
John P. Gyekenyesi of Glenn Research Cen- 
ter, Osama M. Jadaan and Tamas Palfi of 
Ohio Aerospace Institute, Lynn Powers of Case 
Western Reserve University, Stefan Reh of 
ANSYS Inc., and Eric H. Baker of Connecticut 
Reserve Technologies, LLC. Further informa- 
tion is contained in a TSP ( seepage 1 ). 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steve Fedor, Mail 
Stop 4-8, 21 000 Brookpark Road, Cleveland, 
Ohio 44135. Refer to LEW-1 7 682-1/ 4-1. 


$ STRANAL-PMC Version 2.0 

John H. Glenn Research Center, Cleveland, Ohio 


Version 2.0 of the Strain Rate Depen- 
dent Analysis of Polymer Matrix Compos- 
ites (STRANAL-PMC) software has been 
released. A prior version was reported in 
“Analyzing Loads and Strains in Polymer- 
Matrix Composites” (LEW-17227), NASA 
Tech Briefs, Vol. 26, No. 11 (November 
2002), page 36. To recapitulate: Modified 
versions of constitutive equations of vis- 
coplasticity of metals are used to represent 
deformation of a polymeric matrix. The 
equations are applied in a micromechani- 
cal approach, proceeding upward from 
slices of unit cells, through the ply level, to 
the laminate level. The constitutive equa- 


tions are integrated in time by a Runge- 
Kutta technique. To predict the ultimate 
strength of each composite ply, failure cri- 
teria are implemented within the micro- 
mechanics. The inputs to STRANAL-PMC 
are the laminate geometry, properties of 
the fiber and matrix materials, and ap- 
plied stress or strain versus time. The out- 
puts are time-dependent stresses and 
strains at the slice, ply, and laminate levels. 
The improvements in version 2.0 include 
more rigorous representation of hydro- 
static-stress effects in the matrix, refine- 
ment and extension of ply failure models, 
and capabilities to analyze transverse shear 


stresses. Version 2.0 can be implemented 
as a material-model code within transient 
dynamic finite-element codes. 

This program was written by Robert Gold- 
berg and Kelly S. Carney of Glenn Research 
Center, Wieslaw Binienda of the University 
of Akron, andAditi Chattopadhyay of Arizona 
State University. Further information is con- 
tained in a TSP ( seepage 1 ). 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steve Fedor, Mail 
Stop 4-8, 21 000 Brookpark Road, Cleveland, 
Ohio 44135. Refer to LEW-1 7910-1. 


® Micromechanics and Piezo Enhancements of HyperSizer 

John H. Glenn Research Center, Cleveland, Ohio 


The commercial HyperSizer aerospace- 
composite-material-structure-sizing soft- 
ware has been enhanced by incorporating 
capabilities for representing coupled ther- 
mal, piezoelectric, and piezomagnetic ef- 
fects on the levels of plies, laminates, and 
stiffened panels. This enhancement is 
based on a formulation similar to that of 
the pre-existing HyperSizer capability for 
representing thermal effects. As a result of 
this enhancement, the electric and/or 
magnetic response of a material or struc- 
ture to a mechanical or thermal load, or its 
mechanical response to an applied elec- 
tric or magnetic field can be predicted. In 


another major enhancement, a capability 
for representing micromechanical effects 
has been added by establishment of a link- 
age between HyperSizer and Glenn Re- 
search Center’s Micromechanics Analysis 
Code With Generalized Method of Cells 
(MAC/GMC) computer program, which 
was described in several prior NASA Tech 
Briefs articles. The linkage enables Hyper- 
Sizer to localize to the fiber and matrix 
level rather than only to the ply level, mak- 
ing it possible to predict local failures and 
to predict properties of plies from those of 
the component fiber and matrix materi- 
als. Advanced graphical user interfaces 


and database structures have been devel- 
oped to support the new HyperSizer mi- 
cromechanics capabilities. 

These enhancements were made by Steven 
M. Arnold of Glenn Research Center, Brett 
A. Bednarcyk of Ohio Aerospace Institute, and 
Phillip WA. Yarrington and Craig S. Collier 
of Collier Research Corp. Further information 
is contained in a TSP (seepage 1 ). 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed to 
NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steve Fedor, Mail 
Stop 4-8, 21 000 Brookpark Road, Cleveland, 
Ohio 44135. Refer to LEW-17819-1. 
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Single-Phase Rare-Earth Oxide/Aluminum Oxide Glasses 

These glasses are suitable for advanced optical applications. 


Si0 2 and/or Other Oxide(s) 



RE 2 0 3 ► Al 2 0 3 

Mole Percent 


The Colored Area Within the Polygon represents the range of compositions of single-phase glasses de- 
scribed in the text. 


Marshall Space Flight Center, Alabama 

Glasses that comprise rare-earth oxides 
and aluminum oxide plus, optionally, 
lesser amounts of other oxides, have been 
invented. The other oxide (s) can include 
SiO s , B 2 0 3 , Ge0 2 , and/ or any of a variety 
of glass-forming oxides that have been 
used heretofore in making a variety of 
common and specialty glasses. The glasses 
of the invention can be manufactured in 
bulk single-phase forms to ensure near 
uniformity in optical and mechanical 
characteristics, as needed for such devices 
as optical amplifiers, lasers, and optical 
waveguides (including optical fibers). 
These glasses can also be formulated to 
have high indices of refraction, as needed 
in some of such devices. 

The figure presents a ternary phase dia- 
gram showing the range of proportions of 
the main ingredients of several glasses ac- 
cording to the invention. One rare-earth 
oxide that is included in every formulation 
according to the invention is La 2 0 3 , typi- 
cally in an amount between a few mole 
percent and a few tens of mole percent. As 
in the synthesis of other glasses, the oxide 
ingredients are melted together, then the 
melt is cooled sufficiently rapidly to a tem- 
perature below the melting temperatures 
of the crystalline phases and below the liq- 
uidus temperature of the melt. By “suffi- 
ciently rapidly” is meant rapidly enough to 
ensure solidification into glass before crys- 
tallization occurs. The La 2 0 3 helps to pres- 
ent phase separation during solidification. 

Preferably, during solidification, the 
melt should not be allowed to come in con- 
tact with a solid container because such 
contact could give rise to heterogeneous 
nucleation of crystalline material, prevent- 
ing the formation of glass. For making a 
very small amount of glass, contact can be 
prevented by use of a containerless process 
in which the melt is levitated by a gas 
stream and heated by a C0 2 laser, then 
cooled by simply turning off the laser, tn 
another containerless process better suited 
to production of a larger amount of glass, 
the glass melt is levitated on a thin gas film 


formed by gas flowing through a porous 
membrane in a furnace or an electromag- 
netically heated graphite susceptor. 

Alternatively, a relatively large amount of 
glass can be formed in a conventional float 
glass process: First, the oxide ingredients 
are melted together in a crucible. Then the 
glass melt is cast onto a pool of a suitable 
molten metal (e.g., tin or gold) that re- 
mains molten (and, therefore, not crys- 
talline) at a temperature below the temper- 
ature to which the glass melt must be 
rapidly cooled. Another alternative for 
making a relatively large amount of glass is 
to melt the ingredients in a crucible, then 
cast the melt onto a pre-cooled metal die. 


Some heterogeneous nucleation of crystals 
can occur in the glass layer in contact with 
the metal. However, if the rate of cooling is 
fast enough (in effect, if the casting is thin 
enough), then glass is formed before the 
crystals propagate through the thickness of 
the casting. Later, any crystals can be re- 
moved by polishing the glass surface that 
was in contact with the metal die. 

This work was done by J. K. Richard Weber, 
Johan G. Abadie, April D. Hixson, and Paul 
C. Nordine of Containerless Research, Inc., 
for Marshall Space Flight Center. For fur- 
ther information, contact the company at 
cri_info@containerless. com. 

MFS-32 117-1 
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© Tilt/Tip/Piston Manipulator With Base-Mounted Actuators 

The geometry and kinematics of this manipulator would afford advantages for some 
applications. 

Goddard Space Flight Center, Greenbelt, Maryland 



Lengths of Links R-|T 1f R 2 T 2 , and R 3 T 3 are varied to adjust the piston, tilt, and tip coordinates of the platform. 


A proposed three-degree-of-freedom 
(tilt/dp/piston) manipulator, suitable 
for aligning an optical or mechanical 
component, would offer several advan- 
tages over prior such manipulators: 

• Unlike in some other manipulators, no 
actuator would support the weight of 
another actuator: All of the actuators 
would be mounted on a base. Hence, 
there would be less manipulated 
weight. 

• The basic geometry of the manipula- 
tor would afford mechanical advan- 
tage: that is, actuator motions would 
be larger than the motions they pro- 
duce in the manipulated object. Me- 
chanical advantage inherently in- 
creases the accuracy and resolution of 
manipulation. 

• Unlike in some other manipulators, it 
would not be necessary to route power 
and/or data lines through manipula- 
tor joints. 

The proposed manipulator (see fig- 
ure) would include three prismatic actu- 
ators (TiNx, T 2 N 2 , and T 3 N 3 ) mounted 
on the base and operating in the same 
plane. Examples of suitable prismatic ac- 
tuators include lead-screw mechanisms, 
linear hydraulic motors, piezoelectric 
linear drives, inchworm-movement lin- 
ear stepping motors, and linear flexure 
drives. The actuators would control the 
lengths of links RiT 1? R 2 T 2 , and R 3 T 3 . 

Three spherical joints (P 1? P 2 , and P 3 ) 


would be located at the corners of an equi- 
lateral triangle of side length q on the plat- 
form holding the object to be manipu- 
lated. Three inextensible limbs (R 1 P 1 , 
R 2 P 2 , and R 3 P 3 ) having length rwould con- 
nect the spherical joints on the platform 
to revolute joints (R 1? R 2 , and R 3 ) at the 
ends of the actuator-controlled links RiT 1? 
R 2 T 2 , and R 3 T 3 . By varying the lengths of 
these links, one could control the tilt, tip, 
and piston coordinates of the platform. 
Closed-form equations for direct or for- 


ward kinematics of the manipulator (given 
the lengths of the variable links, find the 
tilt, tip, and piston coordinates) have been 
derived. The equations of inverse kinemat- 
ics (find the variable link lengths needed 
to obtain the desired tilt, tip, and piston 
coordinates) have also been derived. 

This work was done by Farhad Tahmasebi of 
Goddard Space Flight Center. Further in- 
formation is contained in a TSP ( seepage 1 ). 
GSC-14874-1 


© Measurement of Model Noise in a Hard-Wall Wind Tunnel 

Spurious noise is suppressed in processing of digitized microphone outputs. 

Ames Research Center, Moffett Field, California 


Identification, analysis, and control of 
fluid-mechanically-generated sound from 
models of aircraft and automobiles in spe- 
cial low-noise, semi-anechoic wind tunnels 
are an important research endeavor. Such 
studies can also be done in aerodynamic 
wind tunnels that have hard walls if phased 
microphone arrays are used to focus on the 


noise-source regions and reject unwanted 
reflections or background noise. Although 
it may be difficult to simulate the total fly- 
over or drive-by noise in a closed wind tun- 
nel, individual noise sources can be iso- 
lated and analyzed. 

An acoustic and aerodynamic study 
was made of a 7-percent-scale aircraft 


model in a NASA Ames 7-by-10-ft (about 
2-by-3-m) wind tunnel for the purpose of 
identifying and attenuating airframe 
noise sources. Simulated landing, take- 
off, and approach configurations were 
evaluated at Mach 0.26. Using a phased 
microphone array mounted in the ceil- 
ing over the inverted model, various 
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Figure 1. The Array of Microphones was 
mounted above the aircraft model in the test 
section of the wind tunnel. A cloth cover has 
been removed from under the microphones to 
make the model visible in this view. 



Figure 2. A Contour Map of Wing-Slat Noise at a 
frequency of 1 2.6 kHz was computed from meas- 
urements made by the microphone array. The 
color-contour range is 8 dB. 

noise sources in the high-lift system, 
landing gear, fins, and miscellaneous 
other components were located and 
compared for sound level and frequency 
at one flyover location. Numerous noise- 
alleviation devices and modifications of 
the model were evaluated. Simultane- 
ously with acoustic measurements, aero- 
dynamic forces were recorded to docu- 
ment aircraft conditions and any 
performance changes caused by geomet- 
ric modifications. 

Most modern microphone-array sys- 
tems function in the frequency domain 
in the sense that spectra of the micro- 


phone outputs are computed, then oper- 
ations are performed on the matrices of 
microphone-signal cross-spectra. The en- 
tire acoustic field at one station in such a 
system is acquired quickly and interro- 
gated during postprocessing. Beam-form- 
ing algorithms are employed to scan a 
plane near the model surface and locate 
noise sources while rejecting most back- 
ground noise and spurious reflections. In 
the case of the system used in this study, 
previous studies in the wind tunnel have 
identified noise sources up to 19 dB 
below the normal background noise of 
the wind tunnel. Theoretical predictions 
of array performance are used to mini- 
mize the width and the side lobes of the 
beam pattern of the microphone array 
for a given test arrangement. 

To capture flyover noise of the in- 
verted model, a 104-element micro- 
phone array in a 622-mm-diameter clus- 
ter was installed in a 19-mm-thick 
poly(methyl methacrylate) plate in the 
ceiling of the test section of the wind 
tunnel above the aircraft model (see Fig- 
ure 1). The microphones were of the 
condenser type, and their diaphragms 
were mounted flush in the array plate, 
which was recessed 12.7 mm into the 
ceiling and covered by a porous aro- 
matic polyamide cloth (not shown in the 
figure) to minimize boundary-layer 
noise. This design caused the level of 
flow noise to be much less than that of 
flush-mount designs. The drawback of 
this design was that the cloth attenuated 
sound somewhat and created acoustic 
resonances that could grow to several dB 
at a frequency of 10 kHz. 

A correction methodology has been 
developed to account for the signal in- 
terference. The first side lobe of the 
beam pattern was 13.4 dB down from 
the peak response at 8 kHz and at an 
angle of 23° from the normal vector: 
these characteristics made it possible to 
obtain good acoustic signals from the 
model when the model was located at a 
distance of 1.11 m from the array. Data 


® Loci-STREAM Version 0.9 

Marshall Space Flight Center, Alabama 


Loci-STREAM is an evolving compu- 
tational fluid dynamics (CFD) software 
tool for simulating possibly chemically 
reacting, possibly unsteady flows in di- 
verse settings, including rocket en- 
gines, turbomachines, oil refineries, 


etc. Loci-STREAM implements a pres- 
sure-based flow-solving algorithm that 
utilizes unstructured grids. (The bene- 
fit of low memory usage by pressure- 
based algorithms is well recognized by 
experts in the field.) The algorithm is 


were acquired at 12,321 scan points in a 
plane encompassing the model. From 
these data, aerodynamic noise from 
sources as small as 6 mm on the model 
surface could be identified easily. 

The microphone signals were digitized 
at a rate of 153,600 samples per second on 
104 channels simultaneously by use of 
analog-to-digital converter circuits and a 
computer. The resulting maximum 
acoustic frequency was 60 kHz with a 
bandwidth of 300 Hz. The data for fre- 
quencies <2 kHz were found to be of mar- 
ginal utility because the microphone 
beam pattern at those frequencies was too 
wide. The data for frequencies >32 kHz 
were found to be of marginal utility be- 
cause at those frequencies, the sources 
were too weak and the side lobes too 
strong. The frequency limits of 2 and 32 
kHz correspond to limits of 140 and 2,240 
Hz, respectively, on the full-scale aircraft. 

A sound-convection correction was in- 
cluded in the processing of the data so 
that sources appeared to come from the 
model rather than being swept down- 
stream. The acoustic sources were de- 
picted, one frequency at a time, as color 
contours on the scan plane with the 
model outline superimposed, as shown 
in Figure 2. Various integration schemes 
have been developed to compute the 
combined effects on a listener and to 
generate narrowband and third-octave 
acoustic spectra. 

Ten airframe noise sources that might 
be important to approach and landing 
noise of the full-scale aircraft were iden- 
tified in the study. The relative strengths 
of these sources and their dependences 
on the configuration of the aircraft were 
documented. Although the data were 
scaled to the frequencies for the full- 
scale aircraft, no extrapolation to full- 
scale flyover was performed. 

This work was done by Paul T. Soderman 
of Ames Research Center. For further in- 
formation, contact the Ames Technology Part- 
nerships Division at (650) 604-2954. 
ARC-14967 


robust for flows at all speeds from zero 
to hypersonic. The flexibility of arbi- 
trary polyhedral grids enables accurate, 
efficient simulation of flows in complex 
geometries, including those of plume- 
impingement problems. The present 
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version — Loci-STREAM version 0.9 — 
includes an interface with the Portable, 
Extensible Toolkit for Scientific Com- 
putation (PETSc) library for access to 
enhanced linear-equation-solving pro- 
grams therein that accelerate conver- 
gence toward a solution. The name 
“Loci” reflects the creation of this soft- 
ware within the Loci computational 


framework, which was developed at 
Mississippi State University for the pri- 
mary purpose of simplifying the writing 
of complex multidisciplinary applica- 
tion programs to run in distributed- 
memory computing environments in- 
cluding clusters of personal computers. 
Loci has been designed to relieve appli- 
cation programmers of the details of 


programming for distributed-memory 
computers. 

This program was written by Jeffrey Wright 
and Siddharth Thakur of Streamline Numer- 
ics, Inc. for Marshall Space Flight Center. 

Further information is contained in a TSP 
(see page 1 ). 

MFS-32303-1 


The Synergistic Engineering Environment 

Langley Research Center, Hampton, Virginia 


The Synergistic Engineering Environ- 
ment (SEE) is a system of software dedi- 
cated to aiding the understanding of 
space mission operations. The SEE can 
integrate disparate sets of data with ana- 
lytical capabilities, geometric models of 
spacecraft, and a visualization environ- 
ment (see figure), all contributing to 
the creation of an interactive simulation 
of spacecraft. Initially designed to satisfy 
needs pertaining to the International 
Space Station, the SEE has been broad- 
ened in scope to include spacecraft 
ranging from those in low orbit around 
the Earth to those on deep-space mis- 
sions. The SEE includes analytical capa- 
bilities in rigid-body dynamics, kinemat- 
ics, orbital mechanics, and payload 
operations. These capabilities enable a 
user to perform real-time interactive en- 
gineering analyses focusing on diverse 
aspects of operations, including flight 
attitudes and maneuvers, docking of vis- 
iting spacecraft, robotic operations, im- 
pingement of spacecraft-engine exhaust 
plumes, obscuration of instrumentation 
fields of view, communications, and al- 
ternative assembly configurations. The 
SEE continues to undergo development 
at Langley Research Center. 



Plumes of Jet Firings can be displayed by the SEE. The user can turn off the firing visualization without 
disabling the visualization of the location. 


This program was written by Jonathan 
Cruz of Langley Research Center and 

Scott Angster of Analytical Mechanics Associ- 


ates, Inc. Further information is contained in 
a TSP ( see page 1 ). 

LAR- 16842-1 


Reconfigurable Software for Controlling Formation Flying 

Goddard Space Flight Center, Greenbelt, Maryland 


Software for a system to control the 
trajectories of multiple spacecraft flying 
in formation is being developed to re- 
flect underlying concepts of (1) a de- 
centralized approach to guidance and 
control and (2) reconfigurability of the 
control system, including reconfigura- 
bility of the software and of control laws. 
The software is organized as a modular 


network of software tasks. The computa- 
tional load for both determining rela- 
tive trajectories and planning maneu- 
vers is shared equally among all 
spacecraft in a cluster. The flexibility 
and robustness of the software are ap- 
parent in the fact that tasks can be 
added, removed, or replaced during 
flight. In a computational simulation of 


a representative formation-flying sce- 
nario, it was demonstrated that the fol- 
lowing are among the services per- 
formed by the software: 

• Uploading of commands from a 
ground station and distribution of the 
commands among the spacecraft, 

• Autonomous initiation and reconfigu- 
ration of formations, 
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• Autonomous formation of teams 
through negotiations among the 
spacecraft, 

• Working out details of high-level 
commands (e.g., shapes and sizes of 
geometrically complex formations) , 


• Implementation of a distributed guidance 
law providing autonomous optimization 
and assignment of target states, and 

• Implementation of a decentralized, 
fuel-optimal, impulsive control law for 
planning maneuvers. 


This work was done by Joseph B. Mueller of 
Princeton Satellite Systems, Inc. for Goddard 
Space Flight Center. Further information 
is contained in a TSP ( see page 1 ). GSC- 
14779-1 


© More About the Tetrahedral Unstructured Software System 

Langley Research Center, Hampton, Virginia 



An Adapted Grid and Flow Solution are shown on the X- 
38 vehicle. This is one of the recent features included in 
TetrUSS. 


Te trUSS is a comprehensive suite of 
computational fluid dynamics (CFD) 
programs that won the Software of the 
Year award in 1996 and has found in- 
creasing use in government, academia, 
and industry for solving realistic flow 
problems (especially in aerodynamics 
and aeroelastics of aircraft having com- 
plex shapes). TetrUSS includes not 
only programs for solving basic equa- 
tions of flow but also programs that af- 
ford capabilities for efficient genera- 
tion and utilization of computational 
grids and for graphical representation 
of computed flows (see figure). The 
2004 version of the Tetrahedral Un- 
structured Software System (TetrUSS), 
which is one of two software systems re- 
ported in “NASA’s 2004 Software of the 
Year,” NASA Tech Briefs, Vol. 28, No. 10 
(October 2004), page 18, has been im- 


proved greatly since 1996. These 
improvements include (1) capa- 
bilities to simulate viscous flow by 
solving the Navier-Stokes equa- 
tions on unstructured grids, (2) 
portability to personal computers 
from diverse manufacturers, (3) 
advanced models of turbulence, 

(4) a parallel-processing version 
of one of the unstructured-grid 
Navier-Stokes-equation-solving 
programs, and (5) advanced pro- 
grams for generating unstruc- 
tured grids. 

These programs were written by 
Khaled S. Abdol-Hamid, Neal T. 

Frink, Craig A. Hunter, Paresh C. 
Parikh, Shaky ar Z. Pizadeh, and Jamshid A. 
Samareh o/'Langley Research Center; Ma- 
haraj K. Bhat of EITI; Mohagna J. Pandya 
of Swales Aerospace; and Matthew J. Grismer 


of the U.S. Air Force Research Laboratory. 
Further information is contained in a TSP 
(see page 1 ). . 

LAR-1 6882-1 


© Computing Flows Using Chimera and Unstructured Grids 

John H. Glenn Research Center, Cleveland, Ohio 

DRAGONFLOW is a computer pro- 
gram that solves the Navier-Stokes equa- 
tions of flows in complexly shaped 
three-dimensional regions discretized 
by use of a direct replacement of arbi- 
trary grid overlapping by nonstructured 
(DRAGON) grid. A DRAGON grid (see 
figure) is a combination of a chimera 
grid (a composite of structured sub- 
grids) and a collection of unstructured 
subgrids. DRAGONFLOW incorpo- 
rates modified versions of two prior 
Navier-Stokes-equation-solving pro- 
grams: OVERFLOW, which is designed 
to solve on chimera grids; and USM3D, 
which is used to solve on unstructured 
grids. A master module controls the in- 
vocation of individual modules in the 
libraries. At each time step of a simu- 
lated flow, DRAGONFLOW is invoked 

on the chimera portion of the This is a DRAGON Grid of the annular turbine cascade. 
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DRAGON grid in alternation with 
USM3D, which is invoked on the un- 
structured subgrids of the DRAGON 
grid. The USM3D and OVERFLOW 
modules then immediately exchange 
their solutions and other data. As a re- 


sult, USM3D and OVERFLOW are cou- 
pled seamlessly. 

This program was written by Meng-Sing 
Liou of Glenn Research Center and Yao 

Zheng of Taitech, Inc. Further information is 
contained in a TSP (see page 1 ). . 


Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steve Fedor, Mail 
Stop 4-8, 21000 Brookpark Road, Cleve- 
land, Ohio 44135. Refer to LEW-17509-1. 


© Avoiding Obstructions in Aiming a High-Gain Antenna 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


The High Gain Antenna Pointing and 
Obstruction Avoidance software per- 
forms computations for pointing a Mars 
Rover high-gain antenna for communi- 
cation with Earth while ( 1 ) avoiding line- 


of-sight obstructions (the Martian terrain 
and other parts of the Rover) that would 
block communication and (2) taking ac- 
count of limits in ranges of motion of an- 
tenna gimbals and of kinematic singular- 


ities in gimbal mechanisms. The software 
uses simplified geometric models of ob- 
structions and of the trajectory of the 
Earth in the Martian sky(see figure). It 
treats all obstructions according to a gen- 
eralized approach, computing and con- 
tinually updating the time remaining be- 
fore interception of each obstruction. In 
cases in which the gimbal-mechanism de- 
sign allows two aiming solutions, the al- 
gorithm chooses the solution that pro- 
vides the longest obstruction-free 
Earth-tracking time. If the communica- 
tion session continues until an obstruc- 
tion is encountered in the current point- 
ing solution and the other solution is 
now unobstructed, then the algorithm 
automatically switches to the other posi- 
tion. This software also notifies commu- 
nication-managing software to cease 
transmission during the switch to the un- 
obstructed position, resuming it when 
the switch is complete. 

This program was written by Khaled Ali 
and Charles Vanelli of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP (seepage 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-42960. 



In this Pointing Strategy Example, the rover is flat and level on the Martian surface. The Hardstop Oc- 
clusions W HA and I/I/hb are indicated as shrouded regions on the celestial sphere. The Earth trajectory 
as a function of time is depicted by e. 


© Analyzing Aeroelastic Stability of a Tilt-Rotor Aircraft 

Langley Research Center, Hampton, Virginia 


Proprotor Aeroelastic Stability Analy- 
sis, now at version 4.5 (PASTA 4.5), is a 
FORTRAN computer program for ana- 
lyzing the aeroelastic stability of a tilt- 
rotor aircraft in the airplane mode of 
flight. The program employs a 10-de- 
gree-of-freedom (DOF), discrete-coordi- 
nate, linear mathematical model of a 
rotor with three or more blades and its 
drive system coupled to a 10-DOF modal 
model of an airframe. The user can se- 
lect which DOFs are included in the 


analysis. Quasi-steady strip-theory aero- 
dynamics is employed for the aerody- 
namic loads on the blades, a quasi-steady 
representation is employed for the aero- 
dynamic loads acting on the vibrational 
modes of the airframe, and a stability-de- 
rivative approach is used for the aerody- 
namics associated with the rigid-body 
DOFs of the airframe. Blade parameters 
that vary with the blade collective pitch 
can be obtained by interpolation from a 
user-defined table. Stability is deter- 


mined by examining the eigenvalues 
that are obtained by solving the coupled 
equations of motions as a matrix eigen- 
value problem. Notwithstanding the rel- 
ative simplicity of its mathematical foun- 
dation, PASTA 4.5 and its predecessors 
have played key roles in a number of en- 
gineering investigations over the years. 

This program was written by Raymond G. 
Kvaternik of Langley Research Center. 
Further information is contained in a TSP 
(see page 1). LAR-17 175-1 
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© Tracking Positions and Attitudes of Mars Rovers 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


The Surface Attitude Position and 
Pointing (SAPP) software, which runs 
on computers aboard the Mars Explo- 
ration Rovers, tracks the positions and 
attitudes of the rovers on the surface of 
Mars. Each rover acquires data on atti- 
tude from a combination of accelerom- 
eter readings and images of the Sun ac- 
quired autonomously, using a 
pointable camera to search the sky for 
the Sun. Depending on the nature of 
movement commanded remotely by 


operators on Earth, the software prop- 
agates attitude and position by use of 
either (1) accelerometer and gyro- 
scope readings or (2) gyroscope read- 
ings and wheel odometry. Where neces- 
sary, visual odometry is performed on 
images to fine-tune the position up- 
dates, particularly on high-wheel-slip 
terrain. The attitude data are used by 
other software and ground-based per- 
sonnel for pointing a high-gain an- 
tenna, planning and execution of driv- 


ing, and positioning and aiming scien- 
tific instruments. 

This work was done by Khaled Ali, Charles 
Vanelli, Jeffrey Biesiadecki, Alejandro San 
Martin, Mark Maimone, Yang Cheng, and 
James Alexander of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP ( see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41876. 


Stochastic Evolutionary Algorithms for Planning Robot Paths 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A computer program implements sto- 
chastic evolutionary algorithms for plan- 
ning and optimizing collision-free paths 
for robots and their jointed limbs. Sto- 
chastic evolutionary algorithms can be 
made to produce acceptably close 
approximations to exact, optimal solu- 
tions for path-planning problems while 
often demanding much less computa- 
tion than do exhaustive-search and de- 
terministic inverse-kinematics algorithms 
that have been used previously for this 
purpose. Hence, the present software is 
better suited for application aboard ro- 
bots having limited computing capabili- 
ties (see figure). The stochastic aspect 
lies in the use of simulated annealing to 
(1) prevent trapping of an optimization 
algorithm in local minima of an energy- 


like error measure by which the fitness of 
a trial solution is evaluated while (2) en- 
suring that the entire multidimensional 



A comparison of Digital Terrain Maps shows reach- 
ability of targets with the FIDO robotic arm. Green 
(light) areas are reachable, with arm path solutions. 
Grey areas are not reachable and red (dark) areas 
indicate no data available for a solution. (Note: 
FIDO is Field Integrated Design and Operations.) 


configuration and parameter space of 
the path-planning problem is sampled 
efficiently with respect to both robot 
joint angles and computation time. Sim- 
ulated annealing is an established tech- 
nique for avoiding local minima in mul- 
tidimensional optimization problems, 
but has not, until now, been applied to 
planning collision-free robot paths by 
use of low-power computers. 

This program was written by Wolfgang Fink, 
Hrand Aghazarian, Terrance Huntsberger, 
and Richard Terrile of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP ( seepage 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-42206. 


© Compressible Flow Toolbox 

John H. Glenn Research Center, Cleveland, Ohio 


The Compressible Flow Toolbox is 
primarily a MATLAB-language imple- 
mentation of a set of algorithms that 
solve approximately 280 linear and 
nonlinear classical equations for com- 
pressible flow. The toolbox is useful for 
analysis of one-dimensional steady flow 
with either constant entropy, friction, 
heat transfer, or Mach number >1. The 
toolbox also contains algorithms for 
comparing and validating the equa- 
tion-solving algorithms against solu- 
tions previously published in open lit- 
erature. The classical equations solved 


by the Compressible Flow Toolbox are 

as follows: 

• The isentropic-flow equations, 

• The Fanno flow equations (pertaining 
to flow of an ideal gas in a pipe with 
friction), 

• The Rayleigh flow equations (pertain- 
ing to frictionless flow of an ideal gas, 
with heat transfer, in a pipe of con- 
stant cross section) , 

• The normal-shock equations, 

• The oblique-shock equations, 

and 

• The expansion equations. 


This program was written by Kevin J. 
Melcher of Glenn Research Center. Fur- 
ther information is contained in a TSP (see 
page 1). 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to: 

NASA Glenn Research Center 
Innovative Partnerships Office 
Attn: Steve Fedor 
Mail Stop 4-8 
21000 Brookpark Road 
Cleveland, Ohio 44135. 

Refer to LEW-17898-1. 
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© Rapid Aeroelastic Analysis of Blade Flutter in Turbomachines 

John H. Glenn Research Center, Cleveland, Ohio 



FIRST MODE, 822 Hz SECOND MODE, 1822 Hz 


Blade Mode Shapes are shown for first and second modes. The first mode is a bending mode at 822 
Hz, and the second mode is a torsion mode at 1,882 Hz. 


The LINFLUX-AE computer code 
predicts flutter and forced responses of 
blades and vanes in turbomachines 
under subsonic, transonic, and super- 
sonic flow conditions. The code solves 
the Euler equations of unsteady flow in a 
blade passage under the assumption 
that the blades vibrate harmonically at 
small amplitudes. The steady-state non- 
linear Euler equations are solved by a 
separate program, then equations for 
unsteady flow components are obtained 
through linearization around the steady- 
state solution. A structural-dynamics 
analysis (see figure) is performed to de- 
termine the frequencies and mode 
shapes of blade vibrations, a preproces- 
sor interpolates mode shapes from the 
structural-dynamics mesh onto the LIN- 
FLUX computational-fluid-dynamics 
mesh, and an interface code is used to 
convert the steady-state flow solution to 
a form required by LINFLUX. Then 
L1NFLUX solves the linearized equa- 
tions in the frequency domain to calcu- 
late the unsteady aerodynamic pressure 
distribution for a given vibration mode, 
frequency, and interblade phase angle. 


A postprocessor uses the unsteady pres- 
sures to calculate generalized aerody- 
namic forces, response amplitudes, and 
eigenvalues (which determine the flut- 
ter frequency and damping) . In compar- 
ison with the TURBO-AE aeroelastic- 
analysis code, which solves the equations 
in the time domain, LINFLUX-AE is 6 to 
7 times faster. 

This program was written by J. J. Trudell, 
O. Mehmed, G. L. Stejko, and M. A. Bakhle 


of Glenn Research Center; T. S. R. Reddy 
of the University of Toledo; M. Montgomery 
of United Technologies; and J. Verdon of Ohio 
Aerospace Institute. Further information is 
contained in a TSP (see page 1 ). 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steve Fedor, Mail 
Stop 4-8, 21000 Brookpark Road, Cleve- 
land, Ohio 44135. Refer to LEW-17880-1. 


© General Flow-Solver Code for Turbomachinery Applications 

Marshall Space Flight Center, Alabama 


Phantom is a computer code intended 
primarily for real-fluid turbomachinery 
problems. It is based on Corsair, an 
ideal-gas turbomachinery code, devel- 
oped by the same authors, which evolved 
from the ROTOR codes from NASA 
Ames. Phantom is applicable to real and 
ideal fluids, both compressible and in- 
compressible, flowing at subsonic, tran- 
sonic, and supersonic speeds. It utilizes 
structured, overset, O- and H-type zonal 
grids to discretize flow fields and repre- 
sent relative motions of components. 


Values on grid boundaries are updated 
at each time step by bilinear interpola- 
tion from adjacent grids. Inviscid fluxes 
are calculated to third-order spatial ac- 
curacy using Roe’s scheme. Viscous 
fluxes are calculated using second-order- 
accurate central differences. The code 
is second-order accurate in time. Turbu- 
lence is represented by a modified Bald- 
win-Lomax algebraic model. The code 
offers two options for determining prop- 
erties of fluids: One is based on equa- 
tions of state, thermodynamic departure 


functions, and corresponding state prin- 
ciples. The other, which is more effi- 
cient, is based on splines generated from 
tables of properties of real fluids. Phan- 
tom currently contains fluid-property 
routines for water, hydrogen, oxygen, ni- 
trogen, kerosene, methane, and carbon 
monoxide as well as ideal gases. 

This work was done by Daniel Dorney of 
Marshall Space Flight Center and Dou- 
glas Sondak of Boston University. Further in- 
formation is contained in a TSP (seepage 1 ). 
MFS-32321-1 


© Code for Multiblock CFD and Heat-Transfer Computations 

John H. Glenn Research Center, Cleveland, Ohio 


The NASA Glenn Research Center 
General Multi-Block Navier-Stokes Con- 
vective Heat Transfer Code, Glenn-HT, 
has been used extensively to predict heat 
transfer and fluid flow for a variety of 


steady gas turbine engine problems. Re- 
cently, the Glenn-HT code has been 
completely rewritten in Fortran 90/95, a 
more object-oriented language that al- 
lows programmers to create code that is 


more modular and makes more efficient 
use of data structures. The new imple- 
mentation takes full advantage of the ca- 
pabilities of the Fortran 90/95 program- 
ming language. As a result, the 
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Glenn-HT code now provides dynamic 
memory allocation, modular design, 
and unsteady flow capability. This allows 
for the heat-transfer analysis of a full tur- 
bine stage. The code has been demon- 
strated for an unsteady inflow condition, 
and gridding efforts have been initiated 
for a full turbine stage unsteady calcula- 
tion. This analysis will be the first to si- 
multaneously include the effects of rota- 
tion, blade interaction, film cooling, and 
tip clearance with recessed tip on tur- 
bine heat transfer and cooling perform- 
ance. Future plans call for the applica- 
tion of the new Glenn-HT code to a 


range of gas turbine engine problems of 
current interest to the heat-transfer 
community. The new unsteady flow ca- 
pability will allow researchers to predict 
the effect of unsteady flow phenomena 
upon the convective heat transfer of tur- 
bine blades and vanes. Work will also 
continue on the development of conju- 
gate heat-transfer capability in the code, 
where simultaneous solution of convec- 
tive and conductive heat-transfer do- 
mains is accomplished. Finally, advanced 
turbulence and fluid flow models and 
automatic gridding techniques are 
being developed that will be applied to 


the Glenn-HT code and solution 
process. 

This program was written by John C. 
Fabian, James D. Heidmann, and Barbara 
L. Lucci of Glenn Research Center; AH A. 
Ameri of the University of Toledo; David L. 
Rigby of QSS, Inc.; and Erlendur Steinthors- 
son A & E Consulting, Inc. Further informa- 
tion is contained in a TSP ( see page 1 ). 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steve Fedor, Mail 
Stop 4-8, 21000 Brookpark Road, Cleve- 
land, Ohio 44135. Refer to LEW-17914-1. 
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Rotating-Pump Design Code 

John H. Glenn Research Center, Cleveland, Ohio 


Pump Design (PUMPDES) is a com- 
puter program for designing a rotating 
pump for liquid hydrogen, liquid oxy- 
gen, liquid nitrogen, water, methane, 
or ethane. Using realistic properties of 
these fluids provided by another pro- 
gram called “GASPAK,” this code per- 
forms a station-by-station, mean-line 
analysis along the pump flow path, ob- 
taining thermodynamic properties of 
the pumped fluid at each station and 
evaluating hydraulic losses along the 
flow path. The variables at each station 
are obtained under constraints that are 
consistent with the underlying physical 


principles. The code evaluates the per- 
formance of each stage and the overall 
pump, fn addition, by judiciously 
choosing the givens and the unknowns, 
the code can perform a geometric in- 
verse design function: that is, it can 
compute a pump geometry that yields a 
closest approximation of given design 
point. The code contains two major 
parts: one for an axial-rotor/inducer 
and one for a multistage centrifugal 
pump. The inducer and the centrifugal 
pump are functionally integrated. The 
code can be used in designing and/or 
evaluating the inducer/centrifugal- 


pump combination or the centrifugal 
pump alone. The code is written in 
standard Fortran 77. 

This program was written by James F. 
Walker and Shu-Cheng Chen of Glenn Re- 
search Center and Dean D. Scheer of Sver- 
drup Technology, Inc. Further information is 
contained in a TSP (see page 1 ). 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steve Fedor, Mail 
Stop 4-8, 21 000 Brookpark Road, Cleveland, 
Ohio 44135. Refer to LEW-17576-1. 
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Manufacturing & Prototyping 


e! Covering a Crucible With Metal Containing Channels 

Metal is deposited on a sacrificial pattern that defines the channels. 

Marshall Space Right Center, Alabama 
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Plastic Tie-Down Tape Wrapped Around Alumina. 
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Wrapped Crucible After Plating With 0.003 Inch of Nickel. 


Part of the Plated Crucible After Removal of the Tape With 
a Piece of Brass Sheet Threaded Through One of the Channels 


Plastic Tie-Down Tape was wrapped around an alumina crucible, the tape-wrapped crucible was 
electroplated with nickel, then the tape was burned out, leaving nickel-walled passages covering 
the crucible. 


In a procedure that partly resembles 
the lost-wax casting process, a crucible 
made of a brittle material (ceramic, 
quartz, or glass) is covered with a layer of 
metal containing channels. The metal 
cover and the channels can serve any or 
all of several purposes, depending upon 
the application: Typically, the metal 
would serve at least partly to reinforce 
the crucible. The channels could be 
used as passages for narrow objects that 
could include thermocouples and heat- 
transfer strips. Alternatively or in addi- 
tion, channels could be used as flow 
paths for liquid or gaseous coolants and 
could be positioned and oriented for po- 
sition- or direction-selective cooling. In 
some cases, the channels could be filled 
with known gases and sealed so that fail- 
ure of the crucibles could be indicated 
by instruments that detect the gases. 

The process consists of three main 
steps. In the first step, a pattern defining 
the channels is formed by wrapping or 
depositing a material in the desired 
channel pattern on the outer surface of 
the crucible. The pattern material can 
be a plastic, wax, low-ash fibrous mate- 
rial, a soluble material, or other suitable 
material that can subsequently be re- 
moved easily. In a proof-of-concept 
demonstration (see figure) , the crucible 
was an alumina cylinder and the mold 
material was plastic tie-down tape. 

In the second step, the patterned cru- 
cible is coated with metal. In one varia- 
tion of the second step, a very thin layer 
containing or consisting of an electri- 
cally conductive material (e.g., gold, sil- 
ver, or carbon) is painted or otherwise 
deposited on the mold-covered cru- 
cible, then the covering metal required 
for the specific application is electrode- 
posited on the very thin conducting 
layer. In another variation of the second 
step, the metal coat is formed by chemi- 
cal vapor deposition. In the proof-of- 
concept demonstration, a layer of nickel 


0.003 in. (-0.08 mm) thick was elec- 
trodeposited. 

In the third step, the patterned mate- 
rial is removed. This is generally done 
by heating the crucible assembly until 
the patterned material melts and runs 
out, vaporizes, and/or decomposes to 
an ash, leaving the channels. Alterna- 
tively, if the patterned material is solu- 
ble, it can be removed by use of a suit- 
able solvent. In the proof-of-concept 
demonstration, the tape was burned 


away by heating the assembly to a tem- 
perature of 600 °C. 

This work was done by Richard N. Grugel of 
Marshall Space Flight Center. Further in- 
formation is contained in a TSP (seepage 1 ). 

This invention has been patented by NASA 
(U.S. Patent No. 6,802,999). Inquiries con- 
cerning nonexclusive or exclusive license for its 
commercial development should be addressed to 
Sammy Nabors, MSFC Commercialization As- 
sistance Lead, at sammy.a.nabors@nasa.gov. 
Refer to MFS-31698-1. 
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Bio-Medical 


© Repairing Fractured Bones by Use of Bioabsorbable Composites 

Less surgery would be necessary, and full strength would be restored sooner. 

Langley Research Center, Hampton, Virginia 


A proposed method of surgical repair of 
fractured bones would incorporate recent 
and future advances in the art of compos- 
ite materials. The composite materials 
used in this method would be biocompati- 
ble and at least partly bioabsorbable: that 
is, during the healing process following 


surgery, they would be wholly or at least 
partly absorbed into the bones and other 
tissues in which they were implanted. Rela- 
tive to the traditional method, the pro- 
posed method would involve less surgery, 
pose less of a risk of infection, provide for 
better transfer of loads across fracture sites, 


and thereby promote better healing while 
reducing the need for immobilization by 
casts and other external devices. 

One requirement that both the tradi- 
tional and proposed methods must satisfy 
is to fix the multiple segments of a broken 
bone in the correct relative positions. Me- 
chanical fixing techniques used in the tra- 
ditional method include the use of plates 
spanning the fracture site and secured to 
the bone by screws, serving of wire along 
the bone across the fracture site, insertion 
of metallic intramedullary rods through 
the hollow portion of the fractured bone, 
and/ or inserting transverse rods through 
the bone, muscle, and skin to stabilize the 
fractured members. After the bone heals, 
a second surgical operation is needed to 
remove the mechanical fixture (s) . In the 
proposed method, there would be no 
need for a second surgical operation. 

The proposed method is based partly on 
the observation that in the fabrication of a 
structural member, it is generally more effi- 
cient and reliable to use multiple small fas- 
teners to transfer load across a joint than to 
use a single or smaller number of larger fas- 
teners, provided that the stress fields of 
neighboring small fasteners do not overlap 
or interact. Also, multiple smaller fasteners 
are more reliable than are larger and fewer 
fasteners. However, there is a trade-off be- 
tween structural efficiency and the cost of 
insertion time and materials. 

The proposed method is further based 
pardy on the conjecture that through-the- 
thickness reinforcements could be excel- 
lent for fixing bone segments for surgical 
repair. The through-the-thickness rein- 
forcements would superficially resemble 
nails in both form and function. Denoted 
small-diameter rods (SDRs) to distinguish 
them from other narrow rods, these rein- 
forcements would be shot or otherwise in- 
serted through adjacent segments of frac- 
tured bone to fix them in their correct 
relative positions (see figure). Shot inser- 
tion would be effected by use an applicator 
that would amount to a miniaturized and 
highly refined version of the pneumatic 
guns often used in carpentry to drive nails 
and brads. The applicator, envisioned to be 
about the size of a ball-point-pen, would be 



BONE WITH COMPOUND FRACTURE, 
FIXED BY USE OF SDRs 



Continuous Wrap Discontinuous Wrap 


BONE WITH COMPOSITE OVERWRAP AFTER 
FIXING WITH USE OF SDRs 


SDRs Would Be Driven Like Nails through adjacent segments of a fractured bone. By use of two or 
more SDRs oriented at different angles, the segments would become locked together. The fractured 
bone would be further immobilized by an overwrapped composite-material shell. 
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driven by pressurized carbon dioxide. To 
further promote stabilization of the seg- 
ments, layers of bone glue could be applied 
to the fracture surfaces prior to insertion 
of the SDRs. The bone glue could be ther- 
apeutically loaded with chemicals to pro- 
mote growth of bone and fight infection. 

SDRs would be produced in a variety of 
diameters to suit specific applications. Typ- 
ical diameters are expected to range from 
0.02 to 0.1 in. (about 0.5 to 2.5 mm). An 
SDR could be fabricated as a matrix/fiber 
composite material: The fibers could be 
made of a biocompatible glass or polymer, 
and the matrix would be made of a bioab- 
sorbable material that could be therapeu- 
tically loaded with bone-growth and infec- 
tion-fighting chemicals. Optionally, the 
fibers could be made of a bioabsorbable 


material, so that in time, nothing of the 
SDR would remain in the healed bone. Yet 
another option would be to fabricate an 
SDR as a metal tube containing bone- 
growth and infection-fighting chemicals 
that would leach out through pores. 

Once the bone segments were fixed in 
place by insertion of SDRs, additional 
structural support and immobilization 
would be provided by wrapping and co- 
curing a multilayer composite-material 
shell around the outer surface of the frac- 
tured bone. This shell would be made of a 
fabric preimpregnated with a curable, 
bioabsorbable matrix resin. From a purely 
structural perspective, it would be prefer- 
able to form the shell around the entire 
circumference of the bone; from a surgi- 
cal perspective, it could be more practical 


to form the shell part way around the cir- 
cumference and strengthen the bone by 
use additional SDRs. 

The curing of the composite would cre- 
ate a rigid structural element integrally 
bonded to the bone. The combination of 
the SDRs and the composite shell may be 
sufficient to restore the bone to its original 
strength, or nearly so. Hence, there would 
be litde or no need for a cast or other ex- 
ternal immobilizing device. At the least, 
the amount of time the patient must wait 
before returning to normal activity should 
be less than it would be if the repair were 
performed by the traditional method. 

This work was done by Gary L. Farley of the U. 
S. Army Research Laboratory for Langley Re- 
search Center. Further information is con- 
tained in a TSP (seepage 1 ). LAR- 16354-1 
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® Kalman Filter for Calibrating a Telescope Focal Plane 

Optimal estimates of scientific and engineering calibration parameters are generated simultaneously. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 



The Focal Plane Is Occupied by two pointing-control reference sensors (PCRSs) and an array of scientific 
instruments. The telescope is maneuvered so that a star is first centroided on PCRS1, then centroided 
on PCRS2, then moved over to the array of scientific instruments to obtain a series of centroid meas- 
urements, then centroided again on PCRS1 . This maneuver makes it possible to attain the main goal of 
the calibration, which is to establish the locations of the instruments in a telescope-pointing reference 
frame. 


The instrument-pointing frame (IPF) 
Kalman filter, and an algorithm that im- 
plements this filter, have been devised for 
calibrating the focal plane of a telescope. 
As used here, “calibration” signifies, more 
specifically, a combination of measure- 
ments and calculations directed toward 
ensuring accuracy in aiming the telescope 
and determining the locations of objects 
imaged in various arrays of photodetec- 
tors in instruments located on the focal 
plane. The IPF Kalman filter was origi- 
nally intended for application to a space- 
borne infrared astronomical telescope, 
but can also be applied to other space- 
borne and ground-based telescopes. 

In the traditional approach to calibra- 
tion of a telescope, (1) one team of ex- 
perts concentrates on estimating param- 
eters (e.g., pointing alignments and 
gyroscope drifts) that are classified as 
being of primarily an engineering na- 
ture, (2) another team of experts con- 
centrates on estimating calibration pa- 
rameters (e.g., plate scales and optical 
distortions) that are classified as being 
primarily of a scientific nature, and (3) 
the two teams repeatedly exchange data 
in an iterative process in which each 
team refines its estimates with the help 
of the data provided by the other team. 
This iterative process is inefficient and 
uneconomical because it is time-con- 
suming and entails the maintenance of 
two survey teams and the development 
of computer programs specific to the re- 
quirements of each team. Moreover, the- 
oretical analysis reveals that the engi- 
neering/ science iterative approach is 
not optimal in that it does not yield the 
best estimates of focal-plane parameters 
and, depending on the application, may 
not even enable convergence toward a 
set of estimates. 

In contrast, in the IPF Kalman-filter 
approach, no attempt is made to distin- 
guish between engineering and scien- 
tific parameters. Hence, there is no 
need for separate engineering and sci- 
entific survey teams, separate software, 
and iteration between the teams. In- 
stead, both engineering and scientific 
focal-plane parameters are estimated to- 


gether, using data taken in the same 
focal-plane survey. The main advantage 
is that the IPF Kalman filter offers 
greater efficiency and economy. In addi- 
tion, the estimates generated by the IPF 
Kalman filter are optimal. 

The IPF Kalman filter is a high-order 
square-root iterated linearized Kalman fil- 
ter, which offers robust numerical condi- 
tioning and a capability to obtain high ac- 
curacy. The filter is parameterized for 
calibrating the focal plane and aligning 
the scientific-instrument photodetector ar- 
rays with respect to the telescope bore- 
sight, all to within a specified tolerance (in 
the original intended application, a focal- 


plane radial standard deviation correspon- 
ding to 0.14 arc second in the sky) . To ob- 
tain this level of accuracy, the filter utilizes 
37 states to estimate desired alignments 
while also correcting for systematic errors 
expected to be caused by optical distor- 
tions, the scale factor and misalignment of 
a scanning mirror, thermomechanically 
induced drifts of alignments among tele- 
scope and instrument frames, and gyro- 
scope bias and bias drift in all axes. 

Other salient features of the IPF 
Kalman filter and algorithm include the 
following: 

• The use of polynomial functions of 
time to characterize such time-depend- 
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ent behaviors as gyroscope drifts and 
thermomechanically induced align- 
ment drifts. Because the polynomial 
coefficients are constant, this feature 
makes it possible to accommodate 
such behaviors while retaining the 
global re-linearization of the Kalman 
filter. 

• A gyroscope-data pre-processing subal- 
gorithm makes it possible to compute 
and store gyroscope sensitivities in ad- 
vance, thereby eliminating the need 
for repeated and time-consuming 
propagation of gyroscope sensitivities 
during each filter cycle. 

• A parameter-masking capability offers 
the option of restricting estimation to 
an arbitrary subset of all the focal- 
plane parameters, thereby affording 


flexibility to match calibration mathe- 
matical models of different levels of fi- 
delity to different scientific instru- 
ments. 

• A “multi-run” feature affords flexibility 
to estimate parameters by use of meas- 
urement data, which has been ac- 
quired during different observing ses- 
sions. 

• An experiment design characterized 
by maneuvers illustrated in the figure 
provides for observability of all desired 
parameters and enables the use of the 
same Kalman filter for a variety of in- 
struments. 

• Measurements using observations of 
both visible and infrared sources can 
be included in the same set of calibra- 
tion data. 


• There is an option to process centroid 
information that is partial in the sense 
that it pertains only to position along a 
single axis of a photodetector array. 
Such information is obtained, for ex- 
ample, when calibrating the entrance 
aperture of spectrometer slit by first 
scanning a source across the narrow 
slit width, and then scanning the 
source along its length. 

• The filter algorithm can be executed 
in one of several optional modes that 
offer compromises between accuracy 
and robustness. 

This work was done by Bryan Kang and 
David Bayard of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP ( see page 1 ). 
NPO-40798 


® Electronic Absolute Cartesian Autocollimator 


Readout is not materially affected by drifts in analog circuitry. 

Goddard Space Flight Center, Greenbelt, Maryland 


An electronic absolute Cartesian auto- 
collimator performs the same basic opti- 
cal function as does a conventional all-op- 
tical or a conventional electronic 
autocollimator but differs in the nature 
of its optical target and the manner in 
which the position of the image of the 
target is measured. The term “absolute” 
in the name of this apparatus reflects the 
nature of the position measurement, 
which, unlike in a conventional elec- 
tronic autocollimator, is based absolutely 
on the position of the image rather than 
on an assumed proportionality between 
the position and the levels of processed 
analog electronic signals. The term 
“Cartesian” in the name of this apparatus 
reflects the nature of its optical target. 

Figure 1 depicts the electronic func- 
tional blocks of an electronic absolute 
Cartesian autocollimator along with its 
basic optical layout, which is the same as 
that of a conventional autocollimator. Re- 
ferring first to the optical layout and func- 
tions only, this or any autocollimator is 
used to measure the compound angular 
deviation of a flat datum mirror with re- 
spect to the optical axis of the autocollima- 
tor itself. The optical components include 
an illuminated target, a beam splitter, an 
objective or collimating lens, and a viewer 
or detector (described in more detail 
below) at a viewing plane. The target and 
the viewing planes are focal planes of the 
lens. Target light reflected by the datum 
mirror is imaged on the viewing plane at 


unit magnification by the collimating lens. 

tf the normal to the datum mirror is 
parallel to the optical axis of the autocol- 
limator, then the target image is centered 
on the viewing plane. Any angular devia- 
tion of the normal from the optical axis 
manifests itself as a lateral displacement 
of the target image from the center. The 
magnitude of the displacement is propor- 
tional to the focal length and to the mag- 
nitude (assumed to be small) of the angu- 
lar deviation. The direction of the 
displacement is perpendicular to the axis 
about which the mirror is slightly tilted. 
Hence, one can determine the amount 
and direction of tilt from the coordinates 
of the target image on the viewing plane. 

tn a conventional all-optical autocolli- 
mator, the target is a first reticle, a tech- 
nician observes the target image 
through an eyepiece, and a second reti- 


cle affixed to the viewing plane is used to 
measure the coordinates of the dis- 
placed image of the first reticle. In a 
conventional electronic autocollimator 
(which could be characterized more ac- 
curately as a conventional optoelec- 
tronic autocollimator), the target is a 
pinhole and a position-sensitive pho- 
todetector is placed at the viewing plane. 
The location of the bright pinhole 
image is measured by use of the posi- 
tion-sensitive photodetector along with 
analog readout circuits. The net outputs 
of these circuits are two sets of voltage 
differences nominally proportional to 
the displacement of the pinhole image 
along two coordinate axes (x and y) in 
the viewing plane. Like all analog de- 
vices and circuits, the position-sensitive 
photodetector and its readout circuits 
exhibit thermal and spontaneous drifts, 



Figure 1. An Electronic Absolute Cartesian Autocollimator includes a conventional autocollimator op- 
tical system with a coded Cartesian-grid target, an image sensor, and a digital image-processing and 
control system. 
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Figure 2. In the Coded Cartesian Grid, each grid cell contains a distinct binary image code that identi- 
fies that cell. The upper code bits in a cell identify the line at the bottom of the cell, while the lower 
bits identify the line at the left of the cell. In the instantaneous field of regard (square box), code bits 
identify vertical lines 2, 3, and 4, and horizontal lines 3, 4, and 5. 


which contribute to errors and lack of 
stability in position measurements. 
Nonuniformity of the position-sensitive 
photodetector also contributes to read- 
out nonlinearity. 

In the electronic absolute Cartesian 
autocollimator, the target is a coded 
Cartesian grid (see Figure 2) and the 
viewing plane is occupied by an image 


sensor. Vertical lines in the target image 
encode azimuthal deflections of the 
datum mirror from the optical axis, 
while horizontal lines encode eleva- 
tional deflections. The planar array of 
pixels of the image sensor intrinsically 
constitutes a fixed high-resolution coor- 
dinate grid. The outputs from the pixels 
are digitized, and the resulting digital 


data are processed to decipher the codes 
in the target image and to determine lo- 
cations of centroids of grid lines, which 
provide angular measurement with a 
granularity nearly one thousand times 
finer than the angular extent of a single 
pixel. Each centroid produces an inde- 
pendent position measurement. Averag- 
ing measurements together naturally in- 
creases readout accuracy and sensitivity. 

The combination of the intrinsic grid 
structure of the image sensor and the 
Cartesian grid of the target image en- 
sures linearity of output and a high de- 
gree of immunity to any non-uniformity 
among responses of individual sensor 
pixels. The coding of the grid ensures 
unambiguous position readout. 

Processing of the target image is not 
subject to drift as a result of weakness of 
signals on the image sensor. At worst, 
weakness of signals increases the pro- 
portion of noise. Therefore, the elec- 
tronic absolute Cartesian autocollima- 
tor includes a servo loop that regulates 
the brightness of illumination to keep 
signal levels optimum. Finally, the elec- 
tronic absolute Cartesian autocollima- 
tor offers one major additional advan- 
tage over a conventional electronic 
autocollimator: The cells of the Carte- 
sian grid effectively constitute a multi- 
plicity of targets that, collectively, makes 
the field of regard of this apparatus 
much larger than that of a conventional 
electronic autocollimator. 

This work was done by Douglas B. Leviton 
of Goddard Space Flight Center. Further 
information is contained in a TSP ( seepage 1 ). 

This invention is owned by NASA, and a 
patent application has been filed. Inquiries 
concerning nonexclusive or exclusive license 
for its commercial development should be ad- 
dressed to the Patent Counsel, Goddard Space 
Flight Center, ( 301 ) 286-7351. Refer to GSC- 
14718-1. 


© Fiber-Optic Gratings for Lidar Measurements of Water Vapor 

These are highly selective, lightweight, tunable optical filters. 

Langley Research Center, Hampton, Virginia 


Narrow-band filters in the form of 
phase-shifted Fabry-Perot Bragg gratings 
incorporated into optical fibers are being 
developed for differential-absorption 
lidar (DIAL) instruments used to meas- 
ure concentrations of atmospheric water 
vapor. The basic idea is to measure the 
relative amounts of pulsed laser light scat- 
tered from the atmosphere at two nearly 


equal wavelengths, one of which coin- 
cides with an absorption spectral peak of 
water molecules and the other correspon- 
ding to no water vapor absorption. As 
part of the DIAL measurement process, 
the scattered light is made to pass 
through a filter on the way to a photode- 
tector. Omitting other details of DIAL for 
the sake of brevity, what is required of the 


filter is to provide a stop band that: 

• Surrounds the water-vapor spectral ab- 
sorption peaks at a wavelength of ~946 
nm, 

• Has a spectral width of at least a couple 
of nanometers, 

• Contains a pass band preferably no 
wider than necessary to accommodate 
the 946.0003-nm-wavelength water- 
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The Transmission Spectrum of a prototype filter was found to contain the desired two narrow pass 
bands within a stop band and was tension-tuned to match the 946-nm water vapor absorption line. 


vapor absorption peak [which has 8.47 
pm full width at half maximum 
(FWHM)], and 

• Contains another pass band at the 
slightly shorter wavelength of 945.9 
nm, where there is scattering of light 
from aerosol particles but no absorp- 
tion by water molecules. 

Whereas filters used heretofore in 
DIAL have had bandwidths of ~300 pm, 
recent progress in the art of fiber-optic 
Bragg-grating filters has made it feasible 
to reduce bandwidths to <20 pm and 
thereby to reduce background noise. An- 
other benefit of substituting fiber-optic 
Bragg-grating filters for those now in use 


would be significant reductions in the 
weights of DIAL instruments. Yet another 
advantage of fiber-optic Bragg-grating fil- 
ters is that their transmission spectra can 
be shifted to longer wavelengths by heat- 
ing or stretching: hence, it is envisioned 
that future DIAL instruments would con- 
tain devices for fine adjustment of trans- 
mission wavelengths through stretching 
or heating of fiber-optic Bragg-grating fil- 
ters nominally designed and fabricated to 
have transmission wavelengths that, in 
the absence of stretching, would be 
slightly too short. 

Prototype fiber-optic Bragg-grating fil- 
ters were designed so that their grating 


structures were chirped and each filter 
included 71-radian phase shifts at two lo- 
cations along its length. In each filter, 
the chirp was characterized by 200 uni- 
form-pitch fields concatenated along a 
total length of about 6 cm. The chirp 
rate was 0.3 nm/cm, with a pitch cen- 
tered at 648.9 nm. The 71-radian phase 
shifts were located at lengthwise posi- 
tions of 29 and 31 cm, respectively. The 
particular combination of chirping pa- 
rameters and phase-shift locations was 
chosen to yield the desired pass bands at 
wavelengths of 945.9 and 946.0003 nm in 
a stop band 2.66 nm wide upon stretch- 
ing of the fiber at a tension equivalent to 
the terrestrial weight of a mass of 140 mg 
(see figure). The filters were fabricated 
in a multistep process, starting with elec- 
tron-beam patterning of step-chirp cor- 
rugations into a mask. Hydrogen-loaded 
single-mode optical fibers were irradi- 
ated through the mask by light from an 
ultraviolet excimer laser, then the fibers 
were annealed by heating. 

The prototype fiber-optic Bragg-grat- 
ing filters were subjected to several tests 
that demonstrated their potential utility 
for DIAL water-vapor measurements. 
Measurements of the transmission spec- 
tra of the filters were found to be well ap- 
proximated by theoretical calculations, 
which were made by use of a piecewise- 
matrix form of a coupled-mode equation. 
Tension tuning was also demonstrated. 

This work was done by Leila B. Vann and 
Russell J. DeYoung of Langley Research 
Center and Stephen J. Mihailov, Ping Lu, 
Dan Grobnic, and Robert Walker of the Com- 
munications Research Centre Canada. Fur- 
ther information is contained in a TSP (see 
pagel). LAR-1 7039-1 


® Simulating Responses of Gravitational-Wave Instrumentation 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Synthetic LISA is a computer pro- 
gram for simulating the responses of the 
instrumentation of the NASA/ESA 
Laser Interferometer Space Antenna 
(LISA) mission, the purpose of which is 
to detect and study gravitational waves. 
Synthetic LISA generates synthetic time 
series of the LISA fundamental noises, 
as filtered through all the time-delay-in- 
terferometry (TDI) observables. (TDI is 
a method of canceling phase noise in 
temporally varying unequal-arm inter- 
ferometers.) Synthetic LISA provides a 
streamlined module to compute the 


TDI responses to gravitational waves, ac- 
cording to a full model of TDI (includ- 
ing the motion of the LISA array and 
the temporal and directional depend- 
ence of the arm lengths). Synthetic 
LISA is written in the C++ programming 
language as a modular package that ac- 
commodates the addition of code for 
specific gravitational wave sources or for 
new noise models. In addition, time se- 
ries for waves and noises can be easily 
loaded from disk storage or electronic 
memory. The package includes a 
Python-language interface for easy, in- 


teractive steering and scripting. 
Through Python, Synthetic LISA can 
read and write data files in Flexible 
Image Transport System (FITS), which 
is a commonly used astronomical data 
format. 

This program was written by John Arm- 
strong, Jeffrey Edlund, and Michele Vallisneri 
of Caltech for NASA’s Jet Propulsion Labo- 
ratory. 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41001. 
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• SOFTC: A Software Correlator for VLBI 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


SOFTC is an advanced software im- 
plementation of a signal correlator for 
very-long-baseline interferometry 
(VLBI) for measuring positions of natu- 
ral celestial objects and distant space- 
craft. Because of increases in processing 
speeds of general-purpose computers, 
software VLBI correlators have become 
viable alternatives to hardware ones. 
The input to SOFTC consists of digi- 
tized samples of raw VLBI-antenna re- 
ceived-signal voltages. Optionally, 
SOFTC also tracks calibration tones su- 
perimposed on the received signals. 


The outputs of SOFTC are (1) phases 
and amplitudes as functions of time and 
frequency for cross-correlated received 
signals and (2) phases and amplitudes 
as functions of time, station, and tone 
number for the calibration tones. 
SOFTC was created to be as accurate as 
possible, capable of processing essen- 
tially any VLBI data, pass strong debug- 
ging tests, have a simple user interface, 
and have no platform dependencies. 
SOFTC is written modularly in the C 
programming language. The great ad- 
vantage of implementing a correlator in 


software, in contradistinction to hard- 
ware, is that it becomes relatively easy 
and much less expensive and time-con- 
suming to adapt, modify, improve, and 
update the correlator. 

This program was written by Stephen Lowe 
of Caltech for NASA’s Jet Propulsion Lab- 
oratory. Further information is contained in 
a TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41072. 


© Progress in Computational Simulation of Earthquakes 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


GeoFEST(P) is a computer program 
written for use in the QuakeSim project, 
which is devoted to development and im- 
provement of means of computational 
simulation of earthquakes. GeoFEST(P) 
models interacting earthquake fault sys- 
tems from the fault-nucleation to the tec- 
tonic scale. The development of Ge- 
oFEST(P) has involved coupling of two 
programs: GeoFEST and the Pyramid 
Adaptive Mesh Refinement Library. Ge- 
oFEST is a message-passing-interface- 
parallel code that utilizes a finite-ele- 
ment technique to simulate evolution of 
stress, fault slip, and plastic/ elastic defor- 


mation in realistic materials like those of 
faulted regions of the crust of the Earth. 
The products of such simulations are 
synthetic observable time-dependent sur- 
face deformations on time scales from 
days to decades. Pyramid Adaptive Mesh 
Refinement Library is a software library 
that facilitates the generation of compu- 
tational meshes for solving physical prob- 
lems. In an application of GeoFEST (P) , a 
computational grid can be dynamically 
adapted as stress grows on a fault. Simu- 
lations on workstations using a few tens 
of thousands of stress and displacement 
finite elements can now be expanded to 


multiple millions of elements with 
greater than 98-percent scaled efficiency 
on over many hundreds of parallel 
processors (see figure). 

This work was done by Andrea Donnel- 
lan, Jay Parker, Gregory Lyzenga, Michele 
Judd, P. Peggy Li, Charles Norton, Edwin 
Tisdale, and Robert Granat of Caltech for 
NASA’s Jet Propulsion Laboratory. 
Further information is contained in a TSP 
(see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41079. 



Models of Landers, CA, Earthquake Deformation are shown at two resolutions. These images show the accuracy improvement going from 82,000 finite el- 
ements on four processors (left) to 1.4 million finite elements on 64 processors (right). 
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Database of Properties of Meteors 

Marshall Space Right Center, Alabama 


A database of properties of meteors, 
and software that provides access to the 
database, are being developed as a con- 
tribution to continuing efforts to model 
the characteristics of meteors with in- 
creasing accuracy. Such modeling is nec- 
essary for evaluation of the risk of pene- 
tration of spacecraft by meteors. For 
each meteor in the database, the record 
will include an identification, date and 
time, radiant properties, ballistic coeffi- 
cient, radar cross section, size, density, 


and orbital elements. The property of 
primary interest in the present case is 
density, and one of the primary goals in 
this case is to derive densities of meteors 
from their atmospheric decelerations. 
The database and software are expected 
to be valid anywhere in the solar system. 
The database will incorporate new data 
plus results of meteoroid analyses that, 
heretofore, have not been readily avail- 
able to the aerospace community. Taken 
together, the database and software con- 


stitute a model that is expected to pro- 
vide improved estimates of densities and 
to result in improved risk analyses for in- 
terplanetary spacecraft. It is planned to 
distribute the database and software on 
a compact disk. 

This program was written by Rob Suggs of 
Marshall Space Flight Center and Coster 
Anthea of Massachusetts Institute of Technol- 
ogy Lincoln Laboratories. For further infor- 
mation, contact Rob Suggs at 
rob.suggs@nasa.gov. MFS-3 2243-1 


® Computing Spacecraft Solar-Cell Damage by Charged Particles 

Goddard Space Flight Center, Greenbelt, Maryland 


General EQFlux is a computer pro- 
gram that converts the measure of the 
damage done to solar cells in outer space 
by impingement of electrons and protons 
having many different kinetic energies 
into the measure of the damage done by 
an equivalent fluence of electrons, each 
having kinetic energy of 1 MeV. Prior to 
the development of General EQFlux, 
there was no single computer program of- 


fering this capability: For a given type of 
solar cell, it was necessary to either per- 
form the calculations manually or to use 
one of three Fortran programs, each of 
which was applicable to only one type of 
solar cell. The problem in developing 
General EQFlux was to rewrite and com- 
bine the three programs into a single pro- 
gram that could perform the calculations 
for three types of solar cells and run in a 


Windows environment with a Windows 
graphical user interface. In comparison 
with the three prior programs, General 
EQFlux is easier to use. 

This program was written by Edward M. 
Gaddy of Goddard Space Flight Center. 
For further information, contact the God- 
dard Innovative Partnerships Office at 
(301) 286-5810. 

GSC-14791-1 


® Thermal Model of a Current-Carrying Wire in a Vacuum 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A computer program implements a ther- 
mal model of an insulated wire carrying 
electric current and surrounded by a vac- 
uum. The model includes the effects of 
Joule heating, conduction of heat along the 
wire, and radiation of heat from the outer 
surface of the insulation on the wire. The 
model takes account of the temperature 
dependences of the thermal and electrical 
properties of the wire, the emissivity of the 


insulation, and the possibility that not only 
can temperature vary along the wire but, in 
addition, the ends of the wire can be ther- 
mally grounded at different temperatures. 
The resulting second-order differential 
equation for the steady-state temperature as 
a function of position along the wire is 
highly nonlinear. The wire is discretized 
along its length, and the equation is solved 
numerically by use of an iterative algorithm 


that utilizes a multidimensional version of 
the Newton-Raphson method. 

This program was written by James Borders 
of Caltech for NASA’s Jet Propulsion Lab- 
oratory. Further information is contained in 
a TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds 
of the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41067. 


® Program for Analyzing Flows in a Complex Network 

Marshall Space Right Center, Alabama 


Generalized Fluid System Simula- 
tion Program (GFSSP) version 4 is a 
general-purpose computer program 
for analyzing steady-state and transient 
flows in a complex fluid network. 
[GFSSP version 2.01 was reported in a 
prior issue of NASA Tech Briefs .] The 


program is capable of modeling com- 
pressibility, fluid transients (e.g., water 
hammers), phase changes, mixtures of 
chemical species, and such externally 
applied body forces as gravitational 
and centrifugal ones. A graphical user 
interface enables the user to interac- 


tively develop a simulation of a fluid 
network consisting of nodes and 
branches. The user can also run the 
simulation and view the results in the 
interface. The system of equations for 
conservation of mass, energy, chemical 
species, and momentum is solved nu- 
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merically by a combination of the New- 
ton-Raphson and successive-substitu- 
tion methods. The program includes 
subroutines that compute thermody- 
namic and thermophysical properties 
for 12 fluids and is integrated with a 
commercial program that gives ther- 
modynamic properties of 36 fluids. 
Eighteen different options are pro- 
vided for modeling momentum 


sources or sinks in the branches. Addi- 
tional capabilities, including new re- 
sistance options, new fluids, and non- 
linear boundary conditions, can be 
added by means of subroutines. An 
audio-visual training CD (compact 
disk) containing lectures, demonstra- 
tion of graphical user interface, and 
tutorial problems is available for learn- 
ing to use the program. 


This program was written by Alok Kumar 
Majumdar of Marshall Space Flight Cen- 
ter. Further information is contained in a 
TSP (see page 1 ). 

This invention is owned by NASA, and a 
patent application has been filed. For further 
information, contact Sammy Nabors, MSFC 
Commercialization Assistance Lead, at 
sammy.a.nabors@nasa.gov. Refer to MFS- 
32125-1. 


® Program Predicts Performance of Optical Parametric Oscillators 

Langley Research Center, Hampton, Virginia 


A computer program predicts the per- 
formances of solid-state lasers that oper- 
ate at wavelengths from ultraviolet 
through mid-infrared and that comprise 
various combinations of stable and unsta- 
ble resonators, optical parametric oscilla- 
tors (OPOs), and sum-frequency genera- 
tors (SFGs), including second-harmonic 
generators (SHGs). The input to the pro- 
gram describes the signal, idler, and 
pump beams; the SFG and OPO crystals; 
and the laser geometry. The program cal- 
culates the electric fields of the idler, 


pump, and output beams at three loca- 
tions (inside the laser resonator, just out- 
side the input mirror, and just outside the 
output mirror) as functions of time for 
the duration of the pump beam. For each 
beam, the electric field is used to calculate 
the fluence at the output mirror, plus 
summary parameters that include the 
centroid location, the radius of curvature 
of the wavefront leaving through the out- 
put mirror, the location and size of the 
beam waist, and a quantity known, vari- 
ously, as a propagation constant or beam- 


quality factor. The program provides a 
typical Windows interface for entering 
data and selecting files. The program can 
include as many as six plot windows, each 
containing four graphs. 

This program was written by Patricia L. 
Cross of Langley Research Center and 
Mark Bowers of Aculight Corporation. Fur- 
ther information is contained in a TSP ( see 
pagel). LAR- 165 29-1 


© Processing TES Level- IB Data 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


TES LIB Subsystem is a computer 
program that performs several func- 
tions for the Tropospheric Emission 
Spectrometer (TES). The term “LIB” 
(an abbreviation of “level IB”), refers to 
data, specific to the TES, on radiometric 
calibrated spectral radiances and their 
corresponding noise equivalent spectral 
radiances (NESRs) , plus ancillary geolo- 
cation, quality, and engineering data. 
The functions performed by TES LIB 
Subsystem include shear analysis, moni- 
toring of signal levels, detection of ice 


build-up, and phase correction and ra- 
diometric and spectral calibration of 
TES target data. Also, the program com- 
putes NESRs for target spectra, writes 
scientific TES level-1 B data to hierarchi- 
cal-data-format (HDF) files for public 
distribution, computes brightness tem- 
peratures, and quantifies interpixel sig- 
nal variability for the purpose of first- 
order cloud and heterogeneous land 
screening by the level-2 software sum- 
marized in the immediately following 
article. This program uses an in-house- 


developed algorithm, called “NUSRT,” 
to correct instrument line-shape factors. 

This program was written by Richard C. De 
Baca, Edwin Sarkissian, Mariyetta Madatyan, 
Douglas Shepard, Scott Gluck, Mark Apolinski, 
James McDuffie, and Dennis Tremblay of Cal- 
tech for N ASA’ s Jet Propulsion Laboratory. 
Further information is contained in a TSP (see 
page 1). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-35218. 


©Automated Camera Calibration 

NASA’s Jet Propulsion Laboratory, Pasadena, California 

Generating calibration data typically in- 
volves taking images of a calibration tar- 
get where the 3D locations of the tar- 
get’s fiducial marks are known, and 
then measuring the 2D locations of the 
fiducial marks in the images. ACAL au- 
tomates the analysis of calibration target 


Automated Camera Calibration 
(ACAL) is a computer program that au- 
tomates the generation of calibration 
data for camera models used in ma- 
chine vision systems. Machine vision 
camera models describe the mapping 
between points in three-dimensional 


(3D) space in front of the camera and 
the corresponding points in two-dimen- 
sional (2D) space in the camera’s 
image. Calibrating a camera model re- 
quires a set of calibration data contain- 
ing known 3D-to-2D point correspon- 
dences for the given camera system. 
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images and greatly speeds the overall 
calibration process. ACAL consists of 
three modules: 

l.ACALDOTS — the primary module 
— takes calibration target images, lo- 
cates and measures the 2D locations 
of the target’s fiducial marks and then 
synthesizes their corresponding 3D lo- 
cations based on knowledge of the cal- 
ibration target’s geometry and its 3D 
location. ACALDOTS handles un- 
even lighting, large-scale variations 
due to range differences, and barrel 
distortion effects of the type found in 
wide-angle lenses. It understands 


both planar and corner-cube (i.e., 
3D) calibration target geometries. 

2. ACALINFO takes the calibration data 
produced by ACALDOTS and esti- 
mates an initial 3D position and orien- 
tation (i.e., camera pose) for the cam- 
era to seed the estimation of the 
remaining camera model parameters. 

3. ACALFIX takes the original calibration 
data and the camera model produced 
from it and generates a refined set of 
calibration data by removing localiza- 
tion errors in the 2D fiducial mark po- 
sitions introduced by perspective fore- 
shortening caused by the calibration 


target’s orientation and geometry. 

The automation in ACAL is robust 
enough that under even extreme image 
conditions, the required input from the 
user averages to no more than about 
one mouse click per target fixture. 

This program was written by Siqi Chen, 
Yang Cheng, and Reg Willson of Caltech for 
NASA’s Jet Propulsion Laboratory. Fur- 
ther information is contained in a TSP (see 
page 1). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41312. 


© Tracking the Martian C0 2 Polar Ice Caps in Infrared Images 

NASA’s Jet Propulsion Laboratory, Pasadena, California 






Dust storms 



Clouds, sublimation 
at the cap edge 


Recession of the Seasonal C0 2 Cap is shown for Northern Martian Spring. 


Researchers at NASA’s Jet Propulsion 
Laboratory have developed a method 
for automatically tracking the polar caps 
on Mars as they advance and recede 
each year (see figure). The seasonal 
Mars polar caps are composed mainly of 
CO s ice and are therefore cold enough 
to stand out clearly in infrared data col- 
lected by the Thermal Emission Imaging 
System (THEMIS) onboard the Mars 
Odyssey spacecraft. The Bimodal Image 
Temperature (BIT) histogram analysis 
algorithm analyzes raw, uncalibrated 
data to identify images that contain both 
“cold” (“polar cap”) and “warm” (“not 
polar cap”) pixels. The algorithm dy- 
namically identifies the temperature 
that separates these two regions. This 
flexibility is critical, because in the ab- 
sence of any calibration, the threshold 
temperature can vary significantly from 
image to image. Using the identified 
threshold, the algorithm classifies each 
pixel in the image as “polar cap” or “not 
polar cap,” then identifies the image row 
that contains the spatial transition from 
“polar cap” to “not polar cap.” While this 
method is useful for analyzing data that 
has already been returned by THEMIS, 
it has even more significance with re- 
spect to data that has not yet been col- 


lected. Instead of seeking the polar cap 
only in specific, targeted images, the 
simplicity and efficiency of this method 
makes it feasible for direct, onboard use. 
That is, THEMIS could continuously 
monitor its observations for any detec- 
tions of the polar-cap edge, producing 
detections over a wide range of spatial 
and temporal conditions. This effort can 
greatly contribute to our understanding 


of long-term climatic change on Mars. 

This work was done by Kiri L. Wagstaff 
Rebecca Castano, and Steve Chien of Caltech 

for NASA’s Jet Propulsion Laboratory. 

Further information is contained in a TSP 
(see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41 732. 


© Processing TES Level-2 Data 

NASA’s Jet Propulsion Laboratory, Pasadena, California 

tics for nadir observations) are also in- 
cluded. The subsystem processes gridded 
meteorological information and extracts 
parameters that can be interpolated to 


TES Level 2 Subsystem is a set of com- 
puter programs that performs functions 
complementary to those of the program 
summarized in the immediately preced- 


ing article. TES Level-2 data pertain to re- 
trieved species (or temperature) profiles, 
and errors thereof. Geolocation, quality, 
and other data (e.g., surface characteris- 
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the appropriate latitude, longitude, and 
pressure level based on the date and 
time. Radiances are simulated using the 
aforementioned meteorological informa- 
tion for initial guesses, and spectroscopic- 
parameter tables are generated. At each 
step of the retrieval, a nonlinear-least- 
squares-solving routine is run over multi- 
ple iterations, retrieving a subset of at- 
mospheric constituents, and error 


analysis is performed. Scientific TES 
Level-2 data products are written in a for- 
mat known as Hierarchical Data Format 
Earth Observing System 5 (HDF-EOS 5) 
for public distribution. 

This software was written by Sassaneh 
Poosti, Sirvard Akopyan, Regina Sakurai, 
Hyejung Yun, Pranjit Saha, Irina Strick- 
land, Kevin Croft, Weldon Smith, Rodney 
Hoffman, John Koffend, Gerard Benenyan, 


9 Smagglce Version 1.8 

John H. Glenn Research Center, Cleveland, Ohio 


Smagglce version 1.8 is a set of soft- 
ware tools for geometrical modeling of, 
and generation of grids that conform to, 
both clean and iced airfoils. A prior ver- 
sion (Smagglce 1.2) was described in 
“Preparing and Analyzing Iced Airfoils” 
(LEW-17399), NASA Tech Briefs, Vol. 28, 
No. 8 (August 2004), page 32. Ice 
shapes, especially those that include 
rough surfaces, pose difficulty in gener- 
ating high-quality grids that are essential 
for predicting airflows by use of compu- 
tational fluid dynamics. Smagglce ver- 
sion 1.8 contains software tools needed 
to overcome this difficulty. For a given 
airfoil, it allows the user to define the 


flow domain, decompose the domain 
into blocks, generate grids, merge grid- 
ded blocks, and control the density and 
smoothness of each grid. Among the 
unique features of version 1.8 is a thin C- 
shaped block, called a “viscous sublayer 
block,” which is wrapped around an iced 
airfoil and its wake line and serves as a 
means to generate highly controlled 
grids near the rough ice surface. Users 
can modify block boundary shapes using 
control points of non-uniform rational 
B-spline (NURBS) curves. Concave ice 
regions can be smoothed during geo- 
metrical modeling or creation of the vis- 
cous sublayer block. 


Hari Nair, Edwin Sarkissian, James Mc- 
Duffie, Ruth Monarrez, David Ho, Benny 
Chan, and Michael Lampel of Caltech for 

NASA’s Jet Propulsion Laboratory. Fur- 
ther information is contained in a TSP ( see 
page 1). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-35212. 


This work was done by Mary B. Vickerman, 
Marivell Baez, Herbert W. Schilling, Barbara 
J. Wilson, Donald C. Braun, Anthony W. 
Hackenberg, James A. Pennline, Rula M. 
Coroneos, and Yung K. Choo of Glenn Re- 
search Center. Further information is con- 
tained in a TSP (see page 1 ). 

Inquiries concerning rights for the com- 
mercial use of this invention should be ad- 
dressed to NASA Glenn Research Center, In- 
novative Partnerships Office, Attn: Steve 
Fedor, Mail Stop 4-8, 21000 Brookpark 
Road, Cleveland, Ohio 44135. Refer to 
LEW-17846-1. 
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Information Sciences 



Solving the Swath Segment Selection Problem 

Several techniques for solving the problem have been tested and compared. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 



TARGET ONLY 





SWATHS WITH EXAMPLES OF SHARDS 


Target Area Is Divided into swaths, segments, and shards for the purpose of mathematical modeling 
in the SSSP. 


Several artificial-intelligence search 
techniques have been tested as means of 
solving the swath segment selection 
problem (SSSP) — a real-world prob- 
lem that is not only of interest in its own 
right, but is also useful as a test bed for 
search techniques in general. In sim- 
plest terms, the SSSP is the problem of 
scheduling the observation times of an 
airborne or spaceborne synthetic-aper- 
ture radar (SAR) system to effect the 
maximum coverage of a specified area 
(denoted the target), given a schedule 
of downlinks (opportunities for radio 
transmission of SAR scan data to a 
ground station), given the limit on the 
quantity of SAR scan data that can be 
stored in an onboard memory between 
downlink opportunities, and given the 
limit on the achievable downlink data 
rate. The SSSP is NP complete (short for 
“nondeterministic polynomial time com- 
plete” — characteristic of a class of in- 
tractable problems that can be solved 
only by use of computers capable of 
making guesses and then checking the 
guesses in polynomial time). 

For the purpose of mathematical 
modeling for the SSSP, a target is ap- 
proximated as being planar and is subdi- 
vided in a pattern of swaths, segments, 
and shards (see figure). A swath is de- 
fined as area, usually rectangular, that 
can be scanned during a given interval 
of time. Segments are portions of swaths 
(also usually rectangular) defined by in- 
tersections of swath and target bound- 
aries. Shards — so named because they 
resemble pieces of broken glass — are 
polygonal areas that are derived from 
segments and that represent pieces of 
target area for which data can be gath- 
ered and downlinked. 

The search techniques tested as 
means of solving the SSSP include the 
following: 

• Forward Dispatch : Segments are added 
in temporal order until the system be- 
comes oversubscribed. When the addi- 
tion of a segment results in an oversub- 
scription, the computation is rolled 
back to remove that segment, then the 
computation resumes without that seg- 


ment. In comparison with the other 
techniques, this technique requires 
the least amount of computation time 
but yields solutions of poor quality. 

• Mixed Integer Programming. Mixed inte- 
ger programming is a hybrid of integer 
and linear programming. It applies to 
a formulation of the SSSP in which the 
relationships among variables and con- 
straints are expressed by matrix-vector 
equations and inequalities. In linear 
programming, the variables can have 
continuous values and one seeks a vec- 
tor of such variables that minimizes or 
maximizes an objective function. An 
integer program includes the addi- 
tional constraint that the variables 
must all have integer values. Integer 
programming gives optimal solutions, 
but requires the greatest amount of 
computation time, even when one ac- 
cepts suboptimal solutions. In mixed 
integer programming, some variables 
are integers and some are continuous. 

• Network Flow Relaxation With Depth-First 


Branch and Bound : One constructs a 
network flow graph — a directed, 
edge-labeled graph that represents the 
flow of information through the prob- 
lem. This graph is used in conjunction 
with depth-first branch-and-bound 
(DFBnB) algorithm. This technique 
provides the same information as does 
integer, mixed integer, or linear pro- 
gramming while requiring less compu- 
tation time. 

The results of tests performed thus far 
have led to the conclusion that for best 
performance, one should start with for- 
ward dispatch followed immediately by 
DFBnB. 

This work was done by Russell Knight and 
Benjamin Smith of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP ( see page 1 ). 

The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPO-40454. 
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© The Spatial Standard Observer 

Degrees of visibility and discriminability of targets in images can be estimated. 

Ames Research Center, Moffett Field, California 



The SSO Was Applied to These Two Images of targets viewed from a UAV. The left image contains a 
target calculated to have a visibility of 4.25 JND; the right image contains a similar target that has 
been reduced in contrast to have a visibility of 1 


The spatial standard observer is a com- 
putational model that provides a measure 
of the visibility of a target in a uniform 
background image or of the visual dis- 
criminability of two images. Standard ob- 
servers have long been used in science 
and industry to quantify the discrim- 
inability of colors. Color standard ob- 
servers address the spectral characteris- 
tics of visual stimuli, while the spatial 
standard observer (SSO) , as its name in- 
dicates, addresses spatial characteristics. 

The SSO is based on a model of human 
vision. The SSO was developed in a 
process that included evaluation of a num- 
ber of earlier mathematical models that 
address optical, physiological, and psy- 
chophysical aspects of spatial characteris- 
tics of human visual perception. Elements 
of the prior models are incorporated into 
the SSO, which is formulated as a compro- 
mise between accuracy and simplicity. The 
SSO operates on a digitized monochrome 
still image or on a pair of such images. The 
SSO consists of three submodels that oper- 
ate sequentially on the input image (s): 

1 . A contrast model, which converts an 
input monochrome image to a lumi- 
nance contrast image, wherein lumi- 
nance values are expressed as excur- 
sions from, and normalized to, a mean; 

2. A contrast-sensitivity-filter model that 
includes an oblique-effect filter (which 
accounts for the decline in contrast sen- 
sitivity at oblique viewing angles) ; and 

3. A spatial summation model, in which 
responses are spatially pooled by rais- 
ing each pixel to the power beta, 
adding the results, and raising the sum 
to the 1/(3 power. In this model, (3=2.9 


was found to be a suitable value. 

The net effect of the SSO is to compute 
a numerical measure of the perceptual 
strength of the single image, or of the vis- 
ible difference (denoted the perceptual 
distance) between two images. The unit 
of a measure used in the SSO is the just 
noticeable difference (JND), which is a 
standard measure of perceptual discrim- 
inability. A target that is just visible has a 
measure of 1 JND. 

The SSO was devised to satisfy an in- 
creasing need for a rapid, objective means 
of estimating degrees of visibility and dis- 
criminability of visual elements in scenes 
observed, not only by humans, but also by 
robotic vision systems, under a variety of 
circumstances. Examples of potential ap- 
plications of the SSO include evaluating 
vision from unpiloted aerial vehicles 
(UAVs) [see figure]; predicting visibility 


JND. 

of UAVs from other aircraft; estimating 
visibility, from a control tower, of aircraft 
on runways; measuring visibility, from a 
distance, of damage on aircraft and on a 
space shutde; evaluation of legibility of 
text, icons, or other symbols; specification 
of resolution of a camera or a display de- 
vice; inspection of display devices during 
manufacturing; estimating the quality of 
compressed digital video imagery; and 
predicting the outcomes of corrective 
laser eye surgery. 

This work was done by Andrew B. Watson 
and Albert J. Ahumada, Jr., of Ames Re- 
search Center. 

This invention is owned by NASA and a 
patent application has been filed. Inquiries 
concerning rights for the commercial use of this 
invention should be addressed to the Ames 
Technology Partnerships Division at (650) 
604-295 4. Refer to ARC-14569-1. 


Less-Complex Method of Classifying MPSK 

Nearly optimal performance can be obtained with less computation. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


An alternative to an optimal method of 
automated classification of signals modu- 
lated with M-ary phase-shift-keying (M-ary 
PSK or MPSK) has been derived. The alter- 
native method is approximate, but it offers 
nearly optimal performance and entails 
much less complexity, which translates to 
much less computation time. 

Modulation classification is becoming 


increasingly important in radio-communi- 
cation systems that utilize multiple data 
modulation schemes and include software- 
defined or software-controlled receivers. 
Such a receiver may “know” little a priori 
about an incoming signal but may be re- 
quired to correctly classify its data rate, 
modulation type, and forward error-cor- 
rection code before properly configuring 


itself to acquire and track the symbol tim- 
ing, carrier frequency, and phase, and ulti- 
mately produce decoded bits. Modulation 
classification has long been an important 
component of military interception of ini- 
tially unknown radio signals transmitted by 
adversaries. Modulation classification may 
also be useful for enabling cellular tele- 
phones to automatically recognize differ- 
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Performances of the ML and Approximate Methods were computed in simulations for a case of non- 
coherent reception of signals that had equal probability of being binary PSK or quaternary PSK. 


ent signal types and configure themselves 
accordingly. 

The concept of modulation classifica- 
tion as outlined in the preceding para- 
graph is quite general. However, at the 
present early stage of development, and 
for the purpose of describing the present 
alternative method, the term “modula- 
tion classification” or simply “classifica- 
tion” signifies, more specifically, a distinc- 
tion between M- ary and M'- ary PSK, 
where M and M' represent two different 
integer multiples of 2. 

Both the prior optimal method and the 
present alternative method require the ac- 
quisition of magnitude and phase values of 
a number (N) of consecutive baseband 
samples of the incoming signal + noise. The 
prior optimal method is based on a maxi- 
mum-likelihood (ML) classification rule 
that requires a calculation of likelihood 
functions for the M and M' hypotheses: 
Each likelihood function is an integral, over 
a full cycle of carrier phase, of a complicated 
sum of functions of the baseband sample 
values, the carrier phase, the carrier-signal 
and noise magnitudes, and M or M\ Then 
the likelihood ratio, defined as the ratio be- 
tween the likelihood functions, is com- 
puted, leading to the choice of whichever 
hypothesis — M or M 1 — is more likely. 

tn the alternative method, the integral 
in each likelihood function is approxi- 
mated by a sum over values of the inte- 


grand sampled at a number, 7, of equally 
spaced values of carrier phase. Used in 
this way, / is a parameter that can be ad- 
justed to trade computational complex- 
ity against the probability of misclassifi- 
cation. tn the limit as / — > oo, one 
obtains the integral form of the likeli- 
hood function and thus recovers the ML 
classification. 

The present approximate method has 
been tested in comparison with the ML 


method by means of computational simu- 
lations. The results of the simulations 
have shown that the performance (as 
quantified by probability of misclassifica- 
tion) of the approximate method is nearly 
indistinguishable from that of the ML 
method (see figure). 

This work was done by Jon Hamkins of Cal- 
tech for NASA’ s Jet Propulsion Laboratory. 
Further information is contained in a TSP (see 
page 1 ). NPO-40965 


Improvement in Recursive Hierarchical Segmentation of Data 

Goddard Space Flight Center, Greenbelt, Maryland 


A further modification has been made 
in the algorithm and implementing soft- 
ware reported in “Modified Recursive Hi- 
erarchical Segmentation of Data” (GSC- 
14681-1), NASA Tech Briefs, Vol. 30, No. 6 
(June 2006), page 51. That software per- 
forms recursive hierarchical segmentation 


of data having spatial characteristics (e.g., 
spectral-image data) . The output of a prior 
version of the software contained artifacts, 
including spurious segmentation-image re- 
gions bounded by processing-window 
edges. The modification for suppressing 
the artifacts, mentioned in the cited article, 



A Segmentation of the Landsat ETM+ Image displayed on the left is shown on the right. The new ap- 
proach eliminates processing artifacts. 


was addition of a subroutine that analyzes 
data in the vicinities of seams to find pairs 
of regions that tend to lie adjacent to each 
other on opposite sides of the seams. 
Within each such pair, pixels in one region 
that are more similar to pixels in the other 
region are reassigned to the other region. 
The present modification provides for a 
parameter ranging from 0 to 1 for control- 
ling the relative priority of merges between 
spatially adjacent and spatially non-adja- 
cent regions. At 1, spatially-adjacent-/spa- 
tially-non-adjacent-region merges have 
equal priority. At 0, only spatially-adjacent- 
region merges (no spectral clustering) are 
allowed. Between 0 and 1, spatially-adja- 
cent-region merges have priority over spa- 
tially-non-adjacent ones. 

This program was written by James C. 
Tilton of Goddard Space Flight Center. 
Further information is contained in a TSP 
(see page 1 ). 

GSC- 14 994-1 
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Using Heaps in Recursive Hierarchical Segmentation of Data 

Goddard Space Flight Center, Greenbelt, Maryland 


A modification to increase processing 
speed has been made in the algorithm 
and implementing software reported in 
“Modified Recursive Hierarchical Seg- 
mentation of Data” (GSC-14681-1), NASA 
Tech Briefs , Vol. 30, No. 6 (June 2006), 
page 51. That software performs recursive 
hierarchical segmentation of data having 
spatial characteristics (e.g., spectral-image 
data) . The segmentation process includes 
an iterative subprocess, in each iteration 
of which it is necessary to determine a 
best pair of regions to merge [merges 


being justified by one or more measure (s) 
similarity of pixels in the regions] . In the 
previously reported version of the algo- 
rithm and software, the choice of a best 
pair of regions to merge involved the use 
of a fully sorted list of regions. That ver- 
sion was computationally inefficient be- 
cause a fully sorted list is not needed: what 
is needed is only the identity of the pair of 
regions characterized by the smallest 
measure of dissimilarity. The present 
modification replaces the use of a fully 
sorted list with the use of data heaps, 


which are computationally more efficient 
for performing the required comparisons 
among dissimilarity measures. The modi- 
fication includes the incorporation of 
standard and modified functions for cre- 
ating and updating data heaps. 

This program was written by James C. 
Tilton of Goddard Space Flight Center. 
Further information is contained in a TSP 
(see page 1 ). 

GSC-14995-1 


Tool for Statistical Analysis and Display of Landing Sites 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


MarsLS is a software tool for analyzing 
statistical dispersion of spacecraft-landing 
sites and displaying the results of its analy- 
ses. Originally intended for the Mars Ex- 
plorer Rover (MER) mission, MarsLS is 
also applicable to landing sites on Earth 
and non-MER sites on Mars. MarsLS is a 
collection of interdependent MATLAB 
scripts that utilize the MATLAB graphical- 
user-interface software environment to 
display landing-site data (see figure) on 
calibrated image-maps of the Martian or 
other terrain. The landing-site data com- 
prise latitude/longitude pairs generated 


by Monte Carlo runs of other computer 
programs that simulate entry, descent, 
and landing. Using these data, MarsLS 
can compute a landing-site ellipse — a 
standard means of depicting the area 
within which the spacecraft can be ex- 
pected to land with a given probability. 
MarsLS incorporates several features for 
the user’s convenience, including capabil- 
ities for drawing lines and ellipses, overlay- 
ing kilometer or latitude/longitude grids, 
drawing and/or specifying lines and/or 
points, entering notes, defining and/or 
displaying polygons to indicate hazards or 


areas of interest, and evaluating haz- 
ardous and/or scientifically interesting 
areas. As part of such an evaluation, 
MarsLS can compute the probability of 
landing in a specified polygonal area. 

This program was written by Geoffrey 
Wawrzyniak, Brian Kennedy, Philip Knocke, 
and John Michel of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP ( see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-35239. 



Martian Landing Site Hazard Maps are shown for Spirit (Gusev) and Opportunity (Meridiani). Here, red = not survivable; green = plains; orange = eroded 
craters; and yellow = subdued craters. 
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Automated Assignment of Proposals to Reviewers 

NASA’s Jet Propulsion Laboratory, Pasadena, California 



The Assignment Process begins with selection of an element from Set 1 and from Set 2 and follows a 
set of rules to avoid biases and other problems in proposal reviews. 


A computer program automates the 
process of selecting unbiased peer review- 
ers of research proposals submitted to 
NASA. Heretofore, such selection has 
been performed by manual searching of 
two large databases subject to a set of as- 
signment rules. One database lists pro- 
posals and proposers; the other database 
lists potential reviewers. The manual 
search takes an average of several weeks 
per proposal. In contrast, the present 
software can perform the selection in sec- 
onds. The program begins by selecting 
one entry from each database, then ap- 
plying the assignment rules to this pair of 
entries. If and only if all the assignment 
rules are satisfied, the chosen reviewer is 
assigned to the chosen proposal. The as- 
signment rules enforced by the program 
are (1) a maximum allowable number of 
proposals assigned to a single reviewer; 
(2) a maximum allowable number of re- 
viewers assigned to a single proposal; (3) 
if the proposing team includes a member 
affiliated with an industry, then the re- 
viewer must not be affiliated with any in- 


dustry; and (4) the reviewer must not be 
a member of the proposing team or affil- 
iated with the same institution as that of a 
member of the proposing team. 

This program was written Faiza Lansing 
and Anil Kantak of Caltech for NASA’s Jet 


Propulsion Laboratory. Further informa- 
tion is contained in a TSP ( see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-40902. 


1 Array-Pattern-Match Compiler for Opportunistic Data Analysis 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A computer program has been written 
to facilitate real-time sifting of scientific 
data as they are acquired to find data pat- 
terns deemed to warrant further analysis. 
The patterns in question are of a type de- 
noted array patterns, which are specified 
by nested parenthetical expressions. 
[One example of an array pattern is 
((>3) 0 (^1)): this pattern matches a vec- 
tor of at least three elements, the first of 
which exceeds 3, the second of which is 
0, and the third of which does not equal 


1.] This program accepts a high-level de- 
scription of a static array pattern and 
compiles a highly optimal and compact 
other program to determine whether 
any given instance of any data array 
matches that pattern. The compiler im- 
plemented by this program is independ- 
ent of the target language, so that as new 
languages are used to write code that 
processes scientific data, they can easily 
be adapted to this compiler. This pro- 
gram runs on a variety of different com- 


puting platforms. It must be run in con- 
junction with any one of a number of 
Lisp compilers that are available com- 
mercially or as shareware. 

This program was written by Mark James 
of Caltech for NASA’s Jet Propulsion Lab- 
oratory. Further information is contained in 
a TSP (see page 1 ). 

This software is available for commercial licens- 
ing. Please contort Karina Edmonds of the Califor- 
nia Institute of Technology at (626) 395-2322. 
Refer to NPO-42096. 


Pre-Processor for Compression of Multispectral Image Data 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A computer program that pre- 
processes multispectral image data has 
been developed to provide the Mars Ex- 
ploration Rover (MER) mission with a 
means of exploiting the additional corre- 
lation present in such data without ap- 
preciably increasing the complexity of 
compressing the data. When used in con- 


junction with ICER, a previously devel- 
oped image-data-compression program, 
this program enables improved compres- 
sion of multispectral images, compared 
to that achievable by use of ICER alone. 
As such, it is a straightforward means of 
achieving much of the gain possible 
from exploiting spectral correlation. 


This preprocessor software accommo- 
dates up to seven images that are differ- 
ent spectral bands of the same scene. 
The software performs an approximate 
discrete cosine transform (DCT) pixel- 
wise across the spectral bands. The soft- 
ware is written for speed; in particular 
the DCT operation performs only inte- 
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ger operations (producing integer out- 
put) and uses multiplications sparingly. 
Separate code is used for each possible 
number of spectral bands, including 
numbers for which fast DCT functions 
are not normally implemented. The 
DCT output is scaled so that, if the origi- 
nal images have a bit depth of at most 12, 
the transformed images are guaranteed 
to have a dynamic range appropriate for 


compression by the ICER software on the 
MER rovers. The resulting transformed 
bands are compressed individually by 
ICER. To reconstruct the images, the 
transformed images are first decom- 
pressed by use of the decompressor for 
ICER, then the resulting reconstructed 
images are passed to an inverse-DCT sub- 
program, which reconstructs the various 
spectral bands. 


This work was done by Matthew Klimesh 
and Aaron Kiely of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP (see page 1 ). 

This software is available for commercial licens- 
ing. Please contact Karina Edmonds of the Cali- 
fornia Institute ofTechnobgy at ( 626) 395-2322. 
Refer to NPO-40835. 


Compressing Image Data While Limiting the Effects of Data Losses 

NASA’s Jet Propulsion Laboratory , Pasadena , California 


ICER is computer software that can per- 
form both lossless and lossy compression 
and decompression of gray-scale-image 
data using discrete wavelet transforms. 
Designed for primary use in transmitting 
scientific image data from distant space- 
craft to Earth, ICER incorporates an 
error-containment scheme that limits the 
adverse effects of loss of data and is well 


suited to the data packets transmitted by 
deep-space probes. The error-contain- 
ment scheme includes utilization of the 
algorithm described in “Partitioning a 
Gridded Rectangle Into Smaller Rectan- 
gles” (NPO-30479), NASA Tech Briefs, Vol. 
28, No. 7 (July 2004), page 56. ICER 
has performed well in onboard compres- 
sion of thousands of images trans- 


mitted from the Mars Exploration Rovers. 

This program was written by Aaron Kiely 
and Matthew Klimesh of Caltech for NASA’s 
Jet Propulsion Laboratory. Further infor- 
mation is contained in a TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-40982. 


Flight Operations Analysis Tool 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Flight Operations Analysis Tool 
(FLOAT) is a computer program that 
partly automates the process of assess- 
ing the benefits of planning spacecraft 
missions to incorporate various combi- 


nations of launch vehicles and pay- 
loads (see figure). Designed primarily 
for use by an experienced systems en- 
gineer, FLOAT makes it possible to 
perform a preliminary analysis of 


trade-offs and costs of a proposed mis- 
sion in days, whereas previously, such 
an analysis typically lasted months. 
FLOAT surveys a variety of prior mis- 
sions by querying data from authorita- 
tive NASA sources pertaining to 20 to 
30 mission and interface parameters 
that define space missions. FLOAT 
provides automated, flexible means 
for comparing the parameters to de- 
termine compatibility or the lack 
thereof among payloads, spacecraft, 
and launch vehicles, and for displaying 
the results of such comparisons. 
Sparseness, typical of the data avail- 
able for analysis, does not confound 
this software. FLOAT effects an itera- 
tive process that identifies modifica- 
tions of parameters that could render 
compatible an otherwise incompatible 
mission set. 

This program was written by Robert Easter, 
Linda Herrell, Richard Pomphrey, James 
Chase, Julie Wertz Chen, Jeffrey Smith, and 
Rebecca Carter of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP ( see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-42306. 



Launch Vehicles 


Spacecraft Buses 


• Payloads 


Secondary spacecraft 
(future capability) 


The Database Overview shows current and future capabilities. 
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Improvement in Visual Target Tracking for a Mobile Robot 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


In an improvement of the visual-tar- 
get-tracking software used aboard a mo- 
bile robot (rover) of the type used to ex- 
plore the Martian surface, an 
affine-matching algorithm has been re- 
placed by a combination of a normal- 
ized-cross-correlation (NCC) algorithm 
and a template-image-magnification al- 
gorithm. Although neither NCC nor 
template-image magnification is new, 
the use of both of them to increase the 
degree of reliability with which features 
can be matched is new. In operation, a 


template image of a target is obtained 
from a previous rover position, then the 
magnification of the template image is 
based on the estimated change in the 
target distance from the previous rover 
position to the current rover position 
(see figure). For this purpose, the target 
distance at the previous rover position is 
determined by stereoscopy, while the 
target distance at the current rover posi- 
tion is calculated from an estimate of the 
current pose of the rover. The template 
image is then magnified by an amount 


corresponding to the estimated target 
distance to obtain a best template image 
to match with the image acquired at the 
current rover position. 

This program was written by Won Kim, 
Adnan Ansar, and Richard Madison of Cal- 
tech for NASA’s Jet Propulsion Labora- 
tory. Further information is contained in a 
TSP ( seepage 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-42682. 



Turn-in-Place Experiments show beginning image (left) and end image (right) after 80° rover rotation. As the rover turns, the mast camera turns in the op- 
posite direction to point to the target. 


Software for Simulating Air Traffic 

Ames Research Center, Moffett Field, California 


Future Air Traffic Management Con- 
cepts Evaluation Tool (FACET) is a sys- 
tem of software for performing compu- 
tational simulations for evaluating 
advanced concepts of advanced air-traf- 
fic management. FACET includes a pro- 
gram that generates a graphical user in- 
terface plus programs and databases 
that implement computational models 
of weather, airspace, airports, navigation 
aids, aircraft performance, and aircraft 
trajectories. Examples of concepts stud- 
ied by use of FACET include aircraft self- 
separation for free flight; prediction of 
air-traffic-controller workload; decision 
support for direct routing; integration 
of spacecraft-launch operations into the 
U.S. national airspace system; and traf- 
fic-flow-management using rerouting, 


metering, and ground delays. Aircraft 
can be modeled as flying along either 
flight-plan routes or great-circle routes 
as they climb, cruise, and descend ac- 
cording to their individual performance 
models. The FACET software is modular 
and is written in the Java and C program- 
ming languages. The architecture of 
FACET strikes a balance between flexi- 
bility and fidelity; as a consequence, 
FACET can be used to model system- 
wide airspace operations over the con- 
tiguous U.S., involving as many as 10,000 
aircraft, all on a single desktop or laptop 
computer running any of a variety of op- 
erating systems. Two notable applica- 
tions of FACET include: (1) reroute con- 
formance monitoring algorithms that 
have been implemented in one of the 


Federal Aviation Administration’s na- 
tionally deployed, real-time, operational 
systems; and (2) the licensing and inte- 
gration of FACET with the commercially 
available Flight Explorer, which is an In- 
ternet-based, real-time flight- tracking 
system. 

This program was written by Banavar 
Sridhar, Karl Bilimoria, and Shon Grabbe of 
Ames Research Center and Gano Chat- 
terji, Kapil Sheth, and Daniel Mulfinger of 
Raytheon Co. Further information is con- 
tained in a TSP (see page 1 ). 

This invention is owned by NASA and a 
patent application has been filed. Inquiries 
concerning rights for the commercial use of 
this invention should be addressed to the 
Ames Technology Partnerships Division at 
(650) 604-2954. Refer to ARC-14653-1. 
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Automated Vectorization of Decision-Based Algorithms 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Virtually all existing vectorization algo- 
rithms are designed to only analyze the 
numeric properties of an algorithm and 
distribute those elements across multiple 
processors. This advances the state of the 
practice because it is the only known sys- 
tem, at the time of this reporting, that 
takes high-level statements and analyzes 
them for their decision properties and 
converts them to a form that allows them 
to automatically be executed in parallel. 


The software takes a high-level source 
program that describes a complex deci- 
sion-based condition and rewrites it as a 
disjunctive set of component Boolean re- 
lations that can then be executed in par- 
allel. This is important because parallel 
architectures are becoming more com- 
monplace in conventional systems and 
they have always been present in NASA 
flight systems. This technology allows 
one to take existing condition-based 


code and automatically vectorize it so it 
naturally decomposes across parallel ar- 
chitectures. 

This program was written by Mark James 
of Caltech for NASA’s Jet Propulsion Lab- 
oratory. Further information is contained in 
a TSP ( see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-42524. 


Grayscale Optical Correlator Workbench 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Grayscale Optical Correlator Work- 
bench (GOCWB) is a computer pro- 
gram for use in automatic target recog- 
nition (ATR). GOCWB performs ATR 
with an accurate simulation of a hard- 
ware grayscale optical correlator (GOC). 
This simulation is performed to test fil- 
ters that are created in GOCWB. Thus, 
GOCWB can be used as a stand-alone 
ATR software tool or in combination 
with GOC hardware for building (target 
training), testing, and optimization of 
filters. The software is divided into three 
main parts, denoted filter, testing, and 


training. The training part is used for as- 
sembling training images as input to a 
filter. The filter part is used for combin- 
ing training images into a filter and opti- 
mizing that filter. The testing part is 
used for testing new filters and for gen- 
eral simulation of GOC output. The cur- 
rent version of GOCWB relies on the 
mathematical software tools from MAT- 
LAB binaries for performing matrix op- 
erations and fast Fourier transforms. 
Optimization of filters is based on an al- 
gorithm, known as OT-MACH, in which 
variables specified by the user are para- 


meterized and the best filter is selected 
on the basis of an average result for cor- 
rect identification of targets in multiple 
test images. 

This program was written by Jay Hanan, 
Hanying Zhou, and Tien-Hsin Chao of Cal- 
tech for NASA’s Jet Propulsion Labora- 
tory. Further information is contained in a 
TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41021. 


© “One-Stop Shopping” for Ocean Remote-Sensing and Model Data 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


OurOcean Portal 2.0 (http:// 
ourocean.jpl.nasa.gov) is a software system 
designed to enable users to easily gain ac- 
cess to ocean observation data, both re- 
mote-sensing and in-situ, configure and 
run an Ocean Model with observation data 
assimilated on a remote computer, and vi- 
sualize both the observation data and the 
model outputs. At present, the observation 
data and models focus on the California 
coastal regions and Prince William Sound 
in Alaska. This system can be used to per- 
form both real-time and retrospective 
analyses of remote-sensing data and model 
outputs. OurOcean Portal 2.0 incorporates 
state-of-the-art information technologies 
(IT) such as MySQL database, Java Web 
Server (Apache/Tomcat) , Live Access 
Server (LAS) , interactive graphics with Java 
Applet at the Client site and MatLab/GMT 
at the server site, and distributed comput- 


ing. OurOcean currendy serves over 20 
real-time or historical ocean data products. 
The data are served in pre-generated plots 
or their native data format. For some of the 
datasets, users can choose different plot- 
ting parameters and produce customized 
graphics. OurOcean also serves 3D Ocean 
Model outputs generated by ROMS (Re- 
gional Ocean Model System) using LAS. 
The Live Access Server (LAS) software, de- 
veloped by the Pacific Marine Environ- 
mental Laboratory (PMEL) of the Na- 
tional Oceanic and Atmospheric 
Administration (NOAA) , is a configurable 
Web-server program designed to provide 
flexible access to geo-referenced scientific 
data. The model output can be views as 
plots in horizontal slices, depth profiles or 
time sequences, or can be downloaded as 
raw data in different data formats, such as 
NetCDF, ASCII, Binary, etc. The interactive 


visualization is provided by graphic soft- 
ware, Ferret, also developed by PMEL. In 
addition, OurOcean allows users with min- 
imal computing resources to configure 
and run an Ocean Model with data assimi- 
lation on a remote computer. Users may se- 
lect the forcing input, the data to be assim- 
ilated, the simulation period, and the 
output variables and submit the model to 
run on a backend parallel computer. 
When the run is complete, the output will 
be added to the LAS server for user to re- 
trieve and examine the results. 

This program was written by P. Peggy Li, 
Quoc Vu, and Yi Chao of Caltech, and Zhi- 
jin Li and Jei-Kook Choi of Raytheon for 
NASA’s Jet Propulsion Laboratory. 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41761. 
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State Analysis Database Tool 

NASA’s Jet Propulsion Laboratory, Pasadena, California 



The State-Based Control Architecture is the foundation of the present software. 


The State Analysis Database Tool soft- 
ware establishes a productive environ- 
ment for collaboration among software 
and system engineers engaged in the de- 
velopment of complex interacting sys- 
tems. The tool embodies State Analysis , a 
model-based system engineering metho- 
dology founded on a state-based control 
architecture (see figure). A state repre- 
sents a momentary condition of an evolv- 
ing system, and a model may describe 
how a state evolves and is affected by 
other states. The State Analysis method- 
ology is a process for capturing system 
and software requirements in the form of 
explicit models and states, and defining 
goal-based operational plans consistent 
with the models. Requirements, models, 
and operational concerns have tradition- 
ally been documented in a variety of sys- 
tem engineering artifacts that address 
different aspects of a mission’s lifecycle. 
In State Analysis, requirements, models, 
and operations information are State 
Analysis artifacts that are consistent and 
stored in a State Analysis Database. The 
tool includes a back-end database, a 
multi-platform front-end client, and 
Web-based administrative functions. The 
tool is structured to prompt an engineer 
to follow the State Analysis methodology, 


to encourage state discovery and model 
description, and to make software re- 
quirements and operations plans consis- 
tent with model descriptions. 

This program was written by Robert Ras- 
mussen and Matthew Bennett of Caltech for 
NASA’s Jet Propulsion Laboratory. Fur- 


ther information is contained in a TSP (see 
page 1). 

The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPO-42869. 


@ Generating CAHV and CAHVOR Images With Shadows in ROAMS 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Part of the Rover Analysis, Modeling 
and Simulation (ROAMS) software 
that synthesizes images of terrain has 
been augmented to make the images 
more realistic. [ROAMS was described 
in “Simulating Operation of a Plane- 
tary Rover” (NPO-30722), NASA Tech 
Briefs, Vol. 28, No. 9 (September 
2004), page 52. ROAMS simulates the 
operation of a robotic vehicle (rover) 
exploring terrain on a remote planet.] 
The images are needed for modeling 
responses of rover cameras that pro- 
vide sensory inputs for machine-vision- 
based algorithms for controlling the 
motion of the rover. The augmented 
image-synthesizing part of the ROAMS 
software supports terrain geometry 
and texture specifiable by the user, 
CAHV and CAHVOR camera models, 
and more-realistic shadowing (see fig- 
ure). (The letters in “CAHV” represent 


vectors in a standard photogrammetric 
model of a pinhole camera. Letters O 
and R in “CAHVOR” represent vectors 
used to model distortions.) A contem- 
plated future version of ROAMS would 
support the CAHVORE model, which 
represents more-general cameras, in- 
cluding those having fish-eye or other 
wide-field-of-view lenses. (Letter E in 
“CAHVORE” represents a vector used 
to model apparent motion of a camera 
entrance pupil.) 

This program was written by Richard 
Madison, Abhinandan Jain, and Marc 
Pomerantz of Caltech for NASA’s Jet Propul- 
sion Laboratory. Further information is 
contained in a TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-42406. 



Examples of Shadowing show terrain and rover 
shadows. Pixels that do not have direct line-of- 
sight to the Sun are darkened. 
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Improving UDP/IP Transmission Without Increasing Congestion 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Datagram Retransmission (DGR) is a 
computer program that, within certain 
limits, ensures the reception of each data- 
gram transmitted under the User Data- 
gram Protocol/Internet Protocol. [User 
Datagram Protocol (UDP) is considered 
unreliable because it does not involve a 
reliability-ensuring connection-initiation 
dialogue between sender and receiver. 
UDP is well suited to issuing of many 
small messages to many different re- 
ceivers.] Unlike prior software for ensur- 
ing reception of UDP datagrams, DGR 
does not contribute to network conges- 
tion by retransmitting data more fre- 


quently as an ever-increasing number of 
messages and acknowledgements is lost. 
Instead, DGR does just the opposite: 
DGR includes an adaptive timeout-inter- 
val-computing component that provides 
maximum opportunity for reception of 
acknowledgements, minimizing retrans- 
mission. By monitoring changes in the 
rate at which message-transmission trans- 
actions are completed, DGR detects 
changes in the level of congestion and re- 
sponds by imposing varying degrees of 
delay on the transmission of new mes- 
sages. In addition, DGR maximizes 
throughput by not waiting for acknowl- 


edgement of a message before sending 
the next message. All DGR communica- 
tion is asynchronous, to maximize effi- 
cient utilization of network connections. 
DGR manages multiple concurrent data- 
gram transmission and acknowledge- 
ment conversations. 

This program was written by Scott Burleigh 
of Caltech for NASA’s Jet Propulsion Labo- 
ratory. Further information is contained in a 
TSP ( seepage 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-40868. 


@ FORTRAN Versions of Reformulated HFGMC Codes 

John H. Glenn Research Center, Cleveland, Ohio 


Several FORTRAN codes have been 
written to implement the reformulated 
version of the high-fidelity generalized 
method of cells (HFGMC). Various as- 
pects of the HFGMC and its predecessors 
were described in several prior NASA 
Tech Briefs articles, the most recent being 
“HFGMC Enhancement of MAC/GMC” 
(LEW-17818-1), NASA Tech Briefs. , Vol. 30, 
No. 3 (March 2006), page 34. The 
HFGMC is a mathematical model of mi- 
cromechanics for simulating stress and 
strain responses of fiber/ matrix and 
other composite materials. The HFGMC 
overcomes a major limitation of a prior 
version of the GMC by accounting for 


coupling of shear and normal stresses 
and thereby affords greater accuracy, al- 
beit at a large computational cost. In the 
reformulation of the HFGMC, the issue 
of computational efficiency was ad- 
dressed: as a result, codes that imple- 
ment the reformulated HFGMC com- 
plete their calculations about 10 times as 
fast as do those that implement the 
HFGMC. The present FORTRAN imple- 
mentations of the reformulated HFGMC 
were written to satisfy a need for compat- 
ibility with other FORTRAN programs 
used to analyze structures and composite 
materials. The FORTRAN implementa- 
tions also afford capabilities, beyond 


those of the basic HFGMC, for modeling 
inelasticity, fiber/matrix debonding, and 
coupled thermal, mechanical, piezo, and 
electromagnetic effects. 

These programs were written by Steven M. 
Arnold of Glenn Research Center and 
Jacob Aboudi and Brett A. Bednarcyk of Ohio 
Aerospace Institute. Further information is 
contained in a TSP ( see page 1 ). . 

Inquiries concerning rights for the com- 
mercial use of this invention should be ad- 
dressed to NASA Glenn Research Center, 
Innovative Partnerships Office, Attn: Steve 
Fedor, Mail Stop 4-8, 21000 Brookpark 
Road, Cleveland, Ohio 44135. Refer to 
LEW-17960-1. 


Program for Editing Spacecraft Command Sequences 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Sequence Translator, Editor, and Ex- 
pander Resource (STEER) is a computer 
program that facilitates construction of se- 
quences and blocks of sequences (here- 
after denoted generally as sequence prod- 
ucts) for commanding a spacecraft. 
STEER also provides mechanisms for 
translating among various sequence prod- 
uct types and quickly expanding activities 
of a given sequence in chronological 
order for review and analysis of the se- 
quence. To date, construction of se- 
quence products has generally been done 
by use of such clumsy mechanisms as text- 
editor programs, translating among se- 


quence product types has been challeng- 
ing, and expanding sequences to time-or- 
dered lists has involved arduous processes 
of converting sequence products to “real” 
sequences and running them through 
Class-A software (defined, loosely, as flight 
and ground software critical to a space- 
craft mission). Also, heretofore, generat- 
ing sequence products in standard for- 
mats has been troublesome because 
precise formatting and syntax are re- 
quired. STEER alleviates these issues by 
providing a graphical user interface con- 
taining intuitive fields in which the user 
can enter the necessary information. The 


STEER expansion function provides a 
“quick and dirty” means of seeing how a 
sequence and sequence block would ex- 
pand into a chronological list, without 
need to use of Class-A software. 

This program was written by Roy Gladden, 
Bruce Waggoner, Mark Kordon, Mahnaz 
Hashemi, David Hanks, and Jose Salcedo of 
Caltech for NASA’s Jet Propulsion Labo- 
ratory. Further information is contained in 
a TSP ( see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41175. 
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© Flight-Tested Prototype of BEAM Software 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Researchers at JPL have completed a 
software prototype of BEAM (Beacon- 
based Exception Analysis for Multi- 
missions) and successfully tested its op- 
eration in flight onboard a NASA 
research aircraft. BEAM (see NASA Tech 
Briefs, Vol. 26, No. 9; and Vol. 27, No. 3) 
is an ISHM (Integrated Systems Health 
Management) technology that automat- 
ically analyzes sensor data and classifies 
system behavior as either nominal or 
anomalous, and further characterizes 
anomalies according to strength, dura- 
tion, and affected signals. BEAM (see 
figure) can be used to monitor a wide 
variety of physical systems and sensor 
types in real time. In this series of tests, 
BEAM monitored the engines of a Dry- 


den Flight Research Center F-18 air- 
craft, and performed onboard, unat- 
tended analysis of 26 engine sensors 
from engine startup to shutdown. The 
BEAM algorithm can detect anomalies 
based solely on the sensor data, which 
includes but is not limited to sensor fail- 
ure, performance degradation, incor- 
rect operation such as unplanned en- 
gine shutdown or flameout in this 
example, and major system faults. 
BEAM was tested on an F-18 simulator, 
static engine tests, and 25 individual 
flights totaling approximately 60 hours 
of flight time. During these tests, BEAM 
successfully identified planned anom- 
alies (in-flight shutdowns of one en- 
gine) as well as minor unplanned anom- 


alies (e.g., transient oil- and fuel-pres- 
sure drops) , with no false alarms or sus- 
pected false-negative results for the pe- 
riod tested. BEAM also detected 
previously unknown behavior in the F- 
18 compressor section during several 
flights. This result, confirmed by direct 
analysis of the raw data, serves as a sig- 
nificant test of BEAM’S capability. 

This program was written by Ryan Mackey, 
Raffi Tikidjian, Mark James, and David 
Wang of Caltech for NASA’s Jet Propulsion 
Laboratory. Further information is con- 
tained in a TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-42727. 



Top-Level BEAM Architecture is used for monitoring physical systems in real time. 


Mission Scenario Development Workbench 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


The Mission Scenario Development 
Workbench (MSDW) is a multidiscipli- 
nary performance analysis software tool 
for planning and optimizing space mis- 
sions. It provides a number of new capa- 
bilities that are particularly useful for 
planning the surface activities on other 
planets. MSDW enables rapid planning 
of a space mission and supports flight- 
system and scientific-instrumentation 
trades. It also provides an estimate of 
the ability of flight, ground, and science 
systems to meet high-level mission goals 
and provides means of evaluating ex- 
pected mission performance at an early 


stage of planning in the project life 
cycle. In MSDW, activity plans and 
equipment-list spreadsheets are inte- 
grated with validated parameterized 
simulation models of spacecraft systems. 
In contrast to traditional approaches in- 
volving worst-case estimates with large 
margins, the approach embodied in 
MSDW affords more flexibility and 
more credible results early in the lifecy- 
cle through the use of validated, vari- 
able-fidelity models of spacecraft sys- 
tems. MSDW is expected to help 
maximize the scientific return on invest- 
ment for space missions by understand- 


ing early the performance required to 
have a successful mission while reducing 
the risk of costly design changes made 
at late stages in the project life cycle. 

This program was designed and written by 
Mark Kordon, John Baker, John Gilbert, 
David Hanks, and Dan Mandutianu of Caltech 
for NASA’s Jet Propulsion Laboratory 
and David Hooper of Emaginit. Further in- 
formation is contained in a TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41382. 
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Marsviewer 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Marsviewer is a multi-platform applica- 
tion designed to aid in quality control, 
browsing, and analysis of original science 
product images (Experiment Data 
Records, or EDRs) and derived image 
data products (Reduced Data Records, or 
RDRs) returned by the Mars Explorer 
Rover (MER) mission. Marsviewer offers 
an abstraction of the products’ organiza- 
tion via a “file finder.” For example, the 
application “understands” the file struc- 
ture and filename conventions of the 
MER Operational Storage Server, help- 
ing the user to navigate this complex file 
system to find desired images. Marsviewer 
also works with a flat file system, remote- 


operations file systems, image-archive file 
systems, and others. All EDRs found for 
a given solar day (Sol) are displayed in a 
list, optionally with thumbnail images. 
Once the user selects an image from the 
list, a tabbed pane conveniently displays 
the original source image and all associ- 
ated RDRs. Marsviewer provides the op- 
tion of overlaying derived images upon 
the source image, resulting in an easier- 
to-interpret color representation of the 
data. Display manipulations such as 
zoom, data range adjustment, contrast 
enhancement, and contour control are 
available. Image metadata (labels) from 
the current image can be displayed and 


searched. The architecture of the pro- 
gram is extensible: new types of RDRs can 
be installed and new file finders can be 
added to adapt the program to different 
file structures and different filename con- 
ventions. This keeps the application flex- 
ible and provides an opportunity for 
reuse with future rover missions. 

This program was written by Nicholas Toole 
and Robert Deen of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-40852. 


Tool for Analysis and Reduction of Scientific Data 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


The Automated Virtual Laboratory Tool 
(AVLT) is designed to be an intelligent sci- 
entific analysis assistant (SAA) system, ded- 
icated to facilitating analysis and reduction 
of data collected by spaceborne scientific 
instruments. Within the SAA, a variety of 
conventional and artificial-intelligence 
software tools are integrated into a uni- 
form system architecture. The AVLT inter- 
faces with the user through a sophisticated 
graphical user interface that is part of the 
SAA environment. Functions of the AVLT 
include the following: 

• Understanding the formats of input 
data files and automatically translating 


the files into whatever other formats are 
required for processing by analysis 
functions provided by users; 

• Providing a graphical workspace, and 
assistance in using the workspace, 
wherein scientists can create computa- 
tional experiments for exploration of 
data, formation of hypotheses, applica- 
tion of analysis functions, interpreta- 
tion, and presentation, using concepts 
that are familiar within the scientists’ 
domains of specialty; 

• Providing a subsystem for planning a 
multistep analysis process to attain a goal 
based partly on prior computational 


experiments or on a priori knowledge; 

• Constructing a knowledge base of 
data-exploration methods and analysis 
and interpretation algorithms; and 

• Providing sophisticated graphical-pres- 
entation software to assist in visualiza- 
tion of data. 

This program was written by Mark James 
of Caltech for NASA’s Jet Propulsion Lab- 
oratory. Further information is contained in 
a TSP (see page 1 ). . 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-42514. 


@ ASPEN Version 3.0 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


The Automated Scheduling and Plan- 
ning Environment (ASPEN) computer 
program has been updated to version 3.0. 
ASPEN as a whole (up to version 2.0) has 
been summarized, and selected aspects of 
ASPEN have been discussed in several pre- 
vious NASA Tech Briefs articles. Restated 
briefly, ASPEN is a modular, reconfig- 
urable, application software framework 
for solving batch problems that involve 
reasoning about time, activities, states, and 
resources. Applications of ASPEN can in- 
clude planning spacecraft missions, sched- 
uling of personnel, and managing supply 


chains, inventories, and production lines. 
ASPEN 3.0 can be customized for a wide 
range of applications and for a variety of 
computing environments that include var- 
ious central processing units and random- 
access memories. Domain-specific reason- 
ing modules (e.g., modules for 
determining orbits for spacecraft) can eas- 
ily be plugged into ASPEN 3.0. Improve- 
ments over other, similar software that 
have been incorporated into ASPEN 3.0 
include a provision for more expressive 
time-line values, new parsing capabilities 
afforded by an ASPEN language based on 


Extensible Markup Language, improved 
search capabilities, and improved inter- 
faces to other, utility-type software (notably 
including MATLAB) . 

This program was written by Gregg Ra- 
bideau, Steve Chien, Russell Knight, Steven 
Schaffer, Daniel Tran, Benjamin Cichy, and 
Robert Sherwood of Caltech for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP (seepage 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41986. 
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@ Secure Display of Space-Exploration Images 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Java EDR Display Interface (JEDI) is 
software for either local display or se- 
cure Internet distribution, to author- 
ized clients, of image data acquired 
from cameras aboard spacecraft en- 
gaged in exploration of remote planets. 
(“EDR” signifies “experimental data 
record,” which, in effect, signifies image 
data.) Processed at NASA’s Multimission 
Image Processing Laboratory (MIPL), 
the data can be from either near-real- 
time processing streams or stored files. 
JEDI uses the Java Advanced Imaging 
application program interface, plus 
input/ output packages that are parts of 


the Video Image Communication and 
Retrieval software of the MIPL, to dis- 
play images. JEDI can be run as either a 
standalone application program or 
within a Web browser as a servlet with an 
applet front end. In either operating 
mode, JEDI communicates using the 
HTTP(s) protocol (s). In the Web- 
browser case, the user must provide a 
password to gain access. For each user 
and/or image data type, there is a con- 
figuration file, called a “personality 
file,” containing parameters that con- 
trol the layout of the displays and the in- 
formation to be included in them. Once 


JEDI has accepted the user’s password, 
it processes the requested EDR (pro- 
vided that user is authorized to receive 
the specific EDR) to create a display ac- 
cording to the user’s personality file. 

This program was written by Cecilia 
Cheng, Gillian Thornhill, and Michael 
McAuley of Caltech for NASA’s Jet Propul- 
sion Laboratory. Further information is 
contained in a TSP ( see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-41246. 
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® Digital Front End for 
Wide-Band VLBI Science 
Receiver 

An upgrade to the very-long-baseline- 
interferometry (VLBI) science receiver 
(VSR) — a radio receiver used in NASA’s 
Deep Space Network (DSN) — is cur- 
rently being implemented. The current 
VSR samples standard DSN intermedi- 
ate-frequency (IF) signals at 256 MHz 
and after digital down-conversion 
records data from up to four 16-MHz 
baseband channels. Currently, IF signals 
are limited to the 265-to-375-MHz range, 
and recording rates are limited to less 
than 80 Mbps. The new digital front 
end, denoted the Wideband VSR, pro- 
vides improvements to enable the re- 
ceiver to process wider bandwidth signals 
and accommodate more data channels 
for recording. The Wideband VSR uti- 
lizes state-of-the-art commercial analog- 
to-digital converter and field-program- 
mable gate array (FPGA) integrated 
circuits, and fiber-optic connections in a 
custom architecture. It accepts IF signals 
from 100 to 600 MHz, sampling the sig- 
nal at 1.28 GHz. The sample data are 
sent to a digital processing module, 
using a fiber-optic link for isolation. The 
digital processing module includes 
boards designed around an Advanced 
Telecom Computing Architecture 
(ATCA) industry-standard backplane. 
Digital signal processing implemented in 
FPGAs down-convert the data signals in 
up to 16 baseband channels with pro- 
grammable bandwidths from 1 kHz to 1 6 
MHz. Baseband samples are transmitted 
to a computer via multiple Ethernet con- 
nections allowing recording to disk at 
rates of up to 1 Gbps. 

This work was done by Andre Jongeling, 
Elliott Sigman, Robert Navarro, Charles 
Goodhart, Steve Rogstad, Kumar Chandra, 
Sue Finley, Joseph Trinh, Melissa Soriano, 
Les White, and Robert Proctor of Caltech and 
Benno Rayhrer (contractor) for NASA’s Jet 
Propulsion Laboratory. Further informa- 
tion is contained in a TSP (see page 1 ). . 
NPO-41 1 91 


Multifunctional Tanks for 
Spacecraft 

A document discusses multifunctional 
tanks as means to integrate additional 
structural and functional efficiencies 


into designs of spacecraft. Whereas 
spacecraft tanks are traditionally de- 
signed primarily to store fluids and only 
secondarily to provide other benefits, 
multifunctional tanks are designed to si- 
multaneously provide multiple primary 
benefits. In addition to one or more 
chamber (s) for storage of fluids, a multi- 
functional tank could provide any or all 
of the following: 

• Passageways for transferring the fluids; 

• Part or all of the primary structure of a 
spacecraft; 

• All or part of an enclosure; 

• Mechanical interfaces to components, 
subsystems, and/ or systems; 

• Paths and surfaces for transferring 
heat; 

• Shielding against space radiation; 

• Shielding against electromagnetic in- 
terference; 

• Electrically conductive paths and sur- 
faces; and 

• Shades and baffles to protect against 
sunlight and/or other undesired light. 
Many different multifunctional-tank 

designs are conceivable. The design of a 
particular tank can be tailored to the re- 
quirements for the spacecraft in which 
the tank is to be installed. For example, 
the walls of the tank can be flat or 
curved or have more complicated 
shapes, and the tank can include an in- 
ternal structure for strengthening the 
tank and/or other uses. 

This work was done by David H. Collins, 
Joseph C. Lewis, and Paul D. MacNeal of 
Caltech for NASA’s Jet Propulsion Labo- 
ratory. Further information is contained in 
a TSP (see page 1 ). 

NPO-41 407 


g! Lightweight, Segmented, 
Mostly Silicon Telescope 
Mirror 

A document presents the concept of a 
curved telescope primary reflector struc- 
ture, made mostly of silicon, that would 
have an areal mass density < 1 kg/ m 2 and 
would be deployed in outer space, where 
it would be operated at a temperature in 
the cryogenic range. The concept pro- 
vides for adjustment of the shape of the 
mirror to maintain the required precise 
optical surface figure despite the flexibil- 
ity inherent in the ultra-lightweight de- 
sign. The structure would include a thin 
mirror layer divided into hexagonal seg- 


ments supported by flexure hinges on a 
lightweight two-layer backing structure. 
Each segment would also be supported at 
three points by sets of piezoelectric linear 
microactuators that could impose small 
displacements along the optical axis. The 
excitations applied to the aforemen- 
tioned microactuators would be chosen 
to effect fine adjustments of the axial po- 
sitions and the orientations of the seg- 
ments relative to the supporting struc- 
ture. Other piezoelectric linear 
microactuators embedded in the backing 
structure would enable control of the dis- 
placements of the segments along the 
hexagonal axes; they would also enable 
control of the curvature of the backing 
structure and, thus, additional control of 
the curvature of the reflector. 

This work was done by Eui-Hyeok Yang of 
Caltech for NASA’s Jet Propulsion Labo- 
ratory. Further information is contained in 
a TSP ( see page 1 ). 

NPO-42106 


€? Assistant for Analyzing 
Tropical-Rain-Mapping 
Radar Data 

A document is defined that describes 
an approach for a Tropical Rain Map- 
ping Radar Data System (TDS). TDS is 
composed of software and hardware ele- 
ments incorporating a two-frequency 
spaceborne radar system for measuring 
tropical precipitation. The TDS would 
be used primarily in generating data 
products for scientific investigations. 
The most novel part of the TDS would 
be expert-system software to aid in the 
selection of algorithms for converting 
raw radar-return data into such primary 
observables as rain rate, path-integrated 
rain rate, and surface backscatter. The 
expert-system approach would address 
the issue that selection of algorithms for 
processing the data requires a signifi- 
cant amount of preprocessing, non-intu- 
itive reasoning, and heuristic applica- 
tion, making it infeasible, in many cases, 
to select the proper algorithm in real 
time. In the TDS, tentative selections 
would be made to enable conversions in 
real time. The expert system would re- 
move straightforwardly convertible data 
from further consideration, and would 
examine ambiguous data, performing 
analysis in depth to determine which al- 
gorithms to select. Conversions per- 
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formed by these algorithms, presumed 
to be correct, would be compared with 
the corresponding real-time conver- 
sions. Incorrect real-time conversions 
would be updated using the correct con- 
versions. 

This work was done by Mark James of Cal- 
tech for NASA’s Jet Propulsion Labora- 
tory. Further information is contained in a 
TSP (see page 1 ). 

The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPO-42515. 


^ Anion-Intercalating 

Cathodes for High-Energy- 
Density Cells 

A report discusses physicochemical is- 
sues affecting a fluoride-intercalating 
cathode that operates in conjunction 


with a lithium ion-intercalating anode in 
a rechargeable electrochemical cell de- 
scribed in a cited prior report. The in- 
stant report also discusses correspond- 
ing innovations made in solvent and 
electrolyte compositions since the prior 
report. The advantages of this cell, rela- 
tive to other lithium-ion-based cells, are 
said to be greater potential (5 V vs. 4 V), 
and greater theoretical cathode specific 
capacity (0.9 to 2.2 A-h/g vs. about 0.18 
A-h/g) . The discussion addresses a need 
for the solvent to be unreactive toward 
the lithium anode and to resist anodic 
oxidation at potentials greater than 
about 4.5 V vs. lithium; the pertinent in- 
novation is the selection of propylene 
carbonate (PC) as a solvent having sig- 
nificantly more stability, relative to other 
solvents that have been tried. The dis- 
cussion also addresses the need for an 
electrolyte additive, denoted an anion 
receptor, to complex the fluoride ion; 


the pertinent innovation is the selection 
of tris(hexafluoroisopropyl) borate as a 
superior alternative to the prior anion 
receptor, which was tris(pentafluo- 
rophenyl) borate. 

This work was done by William West of 
Caltech for NASA’s Jet Propulsion Labo- 
ratory. Further information is contained in 
a TSP (see page 1 ). 

In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Innovative Technology Assets Management 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-4231 6, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 
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